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ABSTRACT 


The  objective  of  this  report  is  to  bring  up  to  date  the  present  state  of  DREV’s 
L(W)WKD  model.  Important  changes  to  the  L(W)WKD  program  have  recently  been  made 
as  a  result  of  exchanges  of  both  models  and  experimental  data  between  Canada  and  France 
during  the  last  few  years.  This  has  led  to  an  agreement  on  a  common  set  of  equations  and 
options  that  are  somewhat  different  from  those  used  in  earlier  versions  of  L(W)WKD  and 
in  earlier  versions  of  the  French  model,  PIRAM. 

The  first  part  of  the  report  discusses  changes  to  the  basic  scaling  equations  that  are 
used  in  the  L(W)WKD  model,  and  the  iterative  method  employed.  It  also  discusses  changes 
to  the  equations  that  relate  the  parameters  produced  by  the  basic  scaling  equations  to 
meteorologically  important  parameters,  and  to  the  model  used  to  calculate  the  refractive 
index  profiles  from  the  visible  to  the  far  infrared.  The  second  part  describes  the  additions 
and  changes  to  the  program’s  modes  of  operation,  its  new  options,  outputs,  and  compilation 
modes.  The  third  part  of  the  report  presents  some  typical  examples  of  the  results  obtainable 
using  the  program  and  a  new  analysis  of  the  MAPTIP  data.  Finally,  some  conclusions  and 
future  improvements  to  the  model  are  discussed. 

RESUME 

Ce  rapport  a  pour  objectif  de  presenter  les  mises  a  jour  recentes  du  module 
L(W)WKD  du  CRDV.  Des  changements  importants  ont  ete  apportes  au  cours  des  demieres 
annees,  principalement  h.  la  suite  du  programme  de  cooperation  entre  le  Canada  et  la  France 
portant  sur  la  modeiisation  de  la  propagation  electromagnetique  a  la  surface  de  la  mer.  En 
particulier,  le  CRDV  et  le  CELAR,  en  France,  se  sont  entendus  sur  une  modelisation 
commune  qui  differe  de  celle  que  Ton  retrouve  dans  la  version  originale  de  leurs  modeles 
respectifs:  L(W)WKD  au  CRDV  et  PIRAM  au  CELAR. 

Dans  la  premiere  partie  du  rapport,  on  presente  les  equations  fondamentales  et  la 
methode  iterative  employee  dans  L(W)WKD  pour  le  calcul  des  parametres  physiques 
intermediaires  qui  decrivent  les  proprietes  de  la  couche  de  surface  et  des  profils  de 
temperature  et  d’humidite  desires.  On  y  discute  aussi  du  calcul  du  profil  d’indice  de 
refractivite  du  visible  jusqu’ a  I’infrarouge  lointain.  Dans  la  deuxieme  partie  du  rapport,  on 
decrit  les  ajouts  et  les  changements  apportes  au  code  informatique:  modes  d’operation  et  de 
compilation,  options  nouvelles,  entrees-sorties,  etc.  Dans  la  demiere  partie  on  presente  des 
exemples  de  calcul  et  une  nouvelle  analyse  des  donnees  de  MAPTIP.  Finalement,  on 
discute  en  conclusion  des  ameliorations  futures  prevues  au  modele. 
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EXECUTIVE  SUMMARY 


Due  to  new  electro-optic  (EO)  sensors,  such  as  IRST  (infrared  search  and  track) 
systems,  that  are  planned  for  use  by  today’s  military,  and  the  need  for  tactical  decision  aids 
(TDAs)  for  such  systems,  the  propagation  of  electromagnetic  radiation  through  the 
atmosphere  is  a  question  of  extreme  importance  and  relevance.  As  the  atmosphere’s 
refractive  index  varies  with  pressure,  temperature  and  water  vapour  content,  all  radiation 
follows  a  curved  path.  Straight  line  propagation  is  only  possible  in  non-refractive 
atmospheres  where  the  index  of  refraction  is  absolutely  uniform.  This  means  that,  in  general, 
the  apparent  elevation  angle  of  a  target  as  seen  by  an  observer  is  not  that  given  by  the  straight 
line  between  the  target  and  the  observer.  This  leads  to  two  phenomena.  The  first  is  the 
phenomenon  commonly  called  sub-refraction.  It  occurs  when  electromagnetic  radiation  from 
a  source  is  bent  away  from  the  surface  of  the  earth.  This  can  lead  to  mirage  effects,  where 
an  observer  can  detect  a  single  source  at  two  different  elevation  angles.  The  second 
phenomenon  is  called  super-refraction  and  occurs  when  electromagnetic  radiation  from  a 
source  is  bent  towards  the  surface  of  the  earth,  such  that  an  observer  can  detect  a  source 
beyond  the  geometrically  defined  horizon.  The  strength  of  these  phenomena  is  dependent 
upon  both  the  meteorological  conditions  and  the  radiation’s  wavelength.  In  the  maritime 
environment,  this  means  that  under  some  conditions  the  distance  to  the  visible  horizon  for 
an  observer  can  be  shorter  than  the  distance  to  the  geometric  horizon  (sub-refraction),  while 
under  other  conditions  it  can  be  further  than  the  geometric  horizon  (super-refraction  or 
ducting).  Knowledge  of  the  distance  to  the  visible  horizon  is  very  important  for  determining 
the  maximum  distance  at  which  an  object  can  be  detected  or  tracked,  and  consequently,  for 
determining  the  maximum  amount  of  time  available  for  making  decisions. 

At  DREV  we  have  been  developing  the  L(W)WKD  model  for  predicting  the  effect  of 
atmospheric  refraction  on  image  distortion,  target  detection  range,  ducting  and  mirage 
formation  in  the  marine  environment  over  the  last  10  years.  Recently,  important  changes 
have  been  made  to  many  of  the  fundamental  equations  that  make  up  the  model  as  a  result  of 
exchanges  between  Canada  and  France.  This  exchange  has  culminated  in  our  agreement  on 
a  common  set  of  equations  and  options  that  have  been  integrated  into  Version  7.09  of 
Canada’s  L(W)WKD  code  and  France’s  PIRAM  code.  Thus,  this  report  is  designed  to 
document  this  commonality  and  to  remark  on  the  changes  from  the  previous  L(W)WKD 
versions.  Changes  to  the  operation  of  the  program  are  also  discussed,  examples  of  the  results 
obtainable  using  the  new  model,  and  a  new  analysis  of  the  MAPTIP  data  are  also  provided 
along  with  some  general  conclusions  and  suggestions  for  future  improvements. 
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1.0  INTRODUCTION 

At  DREV  we  have  been  developing  marine  boundary  layer  (MBL)  models  and 
programs  for  the  last  10  years.  One  such  model,  the  L(Wavy),  Walmsley,  KEL,  DREV  or 
L(W)WKD  model  (Ref.  1),  has  been  developed  to  produce  vertical  profiles  of  the 
refractivity,  the  air  temperature,  the  refractivity  structure  parameter,  and  other  profiles  given 
a  specific  set  of  meteorological  parameters.  These  profiles,  can  then  be  used  by  other 
programs  to  predict  the  imaging  characteristics  of  the  marine  boundary  layer.  For  example, 
refractivity  profiles  can  be  used  by  ray-tracing  programs  to  predict  whether  or  not  an  imager 
will  detect  a  secondary  image  (mirage)  of  a  missile  flying  at  a  height  of  5  m  above  the  wave 
tops  at  a  range  of  15  km  (Ref.  2). 

DGA  (France)  has  also  been  developing  a  MBL  model  called  “Profils  d’Indice  de 
Refraction  en  Atmosphere  Marine”  or  PIRAM  (Ref.  3,  27),  whose  outputs  have  been 
compared  with  those  from  L(W)WKD  since  the  MAPTIP  trial  in  1993  (Ref.  2).  As  a  result 
of  many  technical  exchanges  between  Canada  and  France,  an  agreement  has  recently  been 
obtained  on  a  common  set  of  equations  and  options  that  should  be  used  in  both  the 
L(W)WKD  model  and  the  PIRAM  model. 

This  report  describes  this  common  set  of  equations  and  options,  and  comments  on  the 
changes  this  has  produced  between  the  previous  versions  of  L(W)WKD  (prior  to  Ver.  7)  and 
the  current  version.  The  first  part  of  the  report  discusses  the  basic  scaling  equations  that  are 
derived  from  the  similarity  theory  of  Monin  and  Obukhov  (Ref.  4  &  5)  for  the  near-surface 
atmosphere,  methods  for  their  solution,  and  the  iterative  method  used  to  determine  a 
solution  for  the  L(W)WKD  model.  It  also  presents  the  relationships  between  the 
conservative  variables  used  by  the  basic  scaling  equations  with  the  more  commonly 
measured  meteorological  parameters.  DREV’s  technique  of  incorporating  the  effect  of 
waves  into  the  model  is  also  presented.  Following  this,  it  discusses  the  development  and 
validity  of  a  formula  for  the  calculation  of  the  refractive  index  profiles  from  the  visible  to  the 
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far  infrared  (Refs.  1,  6  &  7).  It  ends  by  presenting  a  method  to  calculate  various  structure 
parameters  that  are  useful  for  determining  the  amount  of  turbulence  within  the  MBL. 

The  second  part  of  the  report  describes,  in  detail,  the  structure  and  operation  of  the 
L(W)WKD  program.  The  program  can  now  be  compiled  into  three  different  executables. 
The  two  BATCH  executables  are  designed  to  be  easily  integrated,  as  separate  modules,  into 
modular  programs  such  as  IRBLEM  (Infrared  Boundary  Layer  Model)(Ref.  8);  while  as  an 
INTERACTIVE  executable,  it  is  designed  to  allow  the  user  to  control  all  aspects  of  the 
program,  and  to  generate  quick  graphical  and  textual  outputs  of  many  parameters  for  further 
analysis.  The  LWWKD  BATCH  executable  produces  vertical  profiles  of  the  refractive 
index,  while  the  WKDCN2  BATCH  executable  produces  vertical  profiles  of  the  refractive 
index  structure  parameter. 

The  report  concludes  with  the  presentation  of  some  typical  examples  obtained  using  the 
L(W)WKD  program,  and  a  new  analysis  of  some  of  the  data  obtained  during  the  Marine 
Aerosol  Properties  and  Thermal  Imager  Performance  (MAPTIP)  trial  (Ref.  2,15) 

Finally,  it  presents  some  conclusions,  and  discusses  various  areas  for  future 
improvements  in  the  L(W)WKD  model. 

This  work  was  carried  out  at  DREV  between  May  1998  and  June  1999  under  work  unit 
lABl  1,  Propagation  Effects  in  the  Marine  Boundary  Layer. 

2.0  LrWiWKD  MARINE  BOUNDARY  LAYER  MODEL 

The  L(W)WKD  model  (Ref.  1)  is  based  on  the  WKD  model  developed  for  DREV  by 
Low  and  Hudak  (Ref.  9)  and  the  Wavy  WKD  (WWKD)  model  developed  at  DREV  by 
Beaulieu  (Ref.  10).  The  objective  in  developing  this  model  was  to  provide  a  simple  and 
effective  technique  for  predicting  the  atmosphere’s  refractivity  profile  in  the  marine  boundary 
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layer  at  several  wavelength  bands  using  a  limited  set  of  meteorological  measurements.  This 
set  includes  a  measurement  of  the  sea  temperature,  and  the  wind  speed  at  a  known  height, 
Zq,  above  the  water  level,  and  measurements  of  the  air  pressure,  air  temperature,  and  relative 
humidity  at  another  known,  but  not  necessarily  different  height,  Zj.  The  model  is  based  on 
the  earlier  work  of  Monin  and  Obukhov  (Refs.  4&  5)  using  similarity  theory  and  on  more 
recent  work  by  Walmsley  (Ref.  11).  The  following  sections  describe  the  model  and  the 
recent  changes  that  have  been  incorporated.  In  this  respect  Appendix  A  is  a  particularly 
important  resource,  as  it  provides  the  reader  with  definitions  for  all  the  defined  constants  and 
variables  used  in  this  report. 


2.1  Scaling  Equations 


As  with  all  MBL  models,  L(W)WKD  assumes  that  the  gradient  of  an  adiabatically 
conservative  quantity,  co,  is  proportional  to  ( l/z)(p^(z/L)  where  z  is  the  height  above  the  water 
level  and  (p^(zfL)  is  a  dimensionless  function.  L(W)WKD  uses  the  three  conservative 
quantities;  wind  speed  (u),  virtual  potential  temperature  (0^),  and  specific  humidity  (q),  with 
the  following  dimensionless  gradient  functions  for  the  wind  speed  ((pu),  virtual  potential 
temperature  ((peJ,  and  specific  humidity  ((p^): 


%(0  = 


KZ  du 

M,  dz 


(la) 

(lb) 

(lc) 


u*,  0y*,  and  q.  are  the  respective  scaling  constants  for  the  wind  speed,  virtual  potential 
temperature  and  specific  humidity,  k  (=  0.4)  is  von  Karman’s  constant,  and  ^  (=  z/L)  is  a 
normalized  height.  L  is  called  the  Monin-Obukhov  length  and  is  defined  by: 


(2) 
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where  g  is  the  acceleration  due  to  gravity,  and  7^^,  (®K)  is  the  average  virtual  air  temperature 
within  the  boundary  layer.  As  required,  L  is  positive  (stable  atmosphere)  when  0^*  is  positive 
and  negative  (unstable  atmosphere)  when  0^.  is  negative.  Furthermore,  L  goes  to  zero  as  u* 
goes  to  zero,  or  as  0^.  goes  to  infinity.  These  formulations  differ  from  that  given  in  Ref.  1 
due  to  the  inclusion  of  von  Karman’s  constant  in  all  three  functions  in  Eq.  1,  and  the  use  of 
the  average  virtual  temperature  rather  than  the  average  temperature  in  Eq.  2.  These 
differences  must  be  kept  in  mind  when  any  comparisons  are  made  between  earlier  versions 
(<  Ver.  7)  of  L(W)WKD  and  later  versions  (^  Ver.  7).  The  functional  form  of  these 
equations  is  now  also  in  agreement  with  that  used  by  the  newest  version  of  PIRAM. 


To  solve  Eq.  1 ,  the  following  gradient  functions  are  generally  used  for  the  wind  speed, 
virtual  potential  temperature,  and  specific  humidity.  The  functions  are  empirically  derived 
from  actual  measurements  and  assume  that  both  the  wind  scaling  constant  (u*)  and  the 
Monin-Obukhov  length  are  non-zero. 


1  +  y,C  /  (1  +  0 

,  c  ^  0  ;  KONDO 

9«(Q  =  - 

1  +  ttjC 

,  C  ^  0  ;  WKD 

(3a) 

(1  -  Pj  0-1/4 

,  C  <  0 

■  1  +  y^C  /  (1  0 

,  C  ^  0  ;  KONDO 

II 

O 

<x> 

9- 

1  +  02^ 

,  C  ^  0  ;  WKD 

(3b) 

(1  -  ur"'" 

,C<o 

1 

'  1  +  Y3C  /  (1  +  0 

,  C  ^  0  ;  KONDO 

II 

o 

s—/ 

9- 

1  +  ttgC 

,  C  ^  0  ;  WKD 

(3c) 

(1  -  PsQ-''" 

,  C  <  0 

where  the  coefficients,  Uj,  Pj,  and  yj  (i  =  1,  2  ,3),  are  constants.  ^  (or  L)  greater  than  zero 
gives  the  appropriate  functions  for  a  stable  atmosphere,  and  ^  (or  L)  less  than  zero  gives  the 
functions  for  an  unstable  atmosphere.  Following  Beaulieu  (Ref.  10),  and  Kondo  (Ref.  12), 
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Eq.  3  is  used  with  the  following  values  for  the  coefficients:  a  =  Uj  =  Oj  =  Uj  =  5, 
P  =  Pi  =  P2  =  Pa  =  16,  and  y  =  y,  =  yj  =  73  =  6.  Two  equations  (KONDO  &  WKD)  are  given 
for  the  stable  case.  The  first  formulation  is  that  originally  proposed  by  Kondo  (Ref.  12)  and 
is  now  the  default  form  used  in  both  PIRAM  and  L(W)WKD.  Earlier  versions  of  L(W)WKD 
used  the  WKD  formulation  proposed  by  KelResearch  (Ref.  9)  as  the  default,  and  it  is  still 
provided  as  an  option  by  L(W)WKD. 


Combining  Eqs.  1  and  3,  and  performing  the  integration  gives  the  following  solutions 
for  u(z),  0y(z)  and  q(z): 

u{z)  =  Mo  +  (M,/K)[ln(z/Zo„)+'P„(z;L)-T^(Zo^;L)]  (4a) 

Mz)  =  Q,o  (0vyK)[ln(z/Zoe )  +  T'g  (z;L)  -  (4b) 

<i(z)  =  qo+ (qA)[Hz/zo^)  +  ^g(r,L)  -  ^^(zo,;^)]  (4c) 


where  the  integrated  gradient  functions,  'F,  are  given  by 


'P„(z;L) 


y  ln(l  +0  ,  C  s  0  ;  KONDO 

aC  ,  C  ^  0  ;  WKD 

2tan"*j:  -  21n(l  +x)  -  ln(l  +x^)  ,  ^  <  0 


(5a) 


'Fe  (z;L)  =  'F^fz;!)  - 


y  ln(l  +  0 
■  aC 
-21n(l  +x^) 


,  c  ^  0  ;  KONDO 
,  C  ^  0  ;  WKD 
.  C  <  0 


(5b) 


and 


X  =  (1-pcr. 


(5c) 


Zqu,  Zoe„  and  Zoq  are  the  respective  roughness  lengths  for  the  wind  speed,  virtual  potential 
temperature  and  specific  humidity.  Uq,  0vo,  and  qo  are  constants  of  integration  for  the  wind 
speed,  virtual  potential  temperature  and  specific  humidity  evaluated  at  their  respective 
roughness  length.  This  set  of  equations  is  composed  of  three  equations  with  10  unknowns 


UNCLASSMED 

6 


(including  L).  Before  presenting  any  methods  for  their  determination,  it  is  important  that  we 
first  define  some  important  meteorological  concepts  and  their  relationships. 


Before  discussing  definitions  of  potential  temperature,  virtual  temperature,  water  vapour 
pressure,  etc.,  it  is  important  to  conunent  on  the  conceptual  reason  for  many  of  the  changes 
that  have  been  made  with  respect  to  those  presented  in  the  previous  L(W)WKD  report 
(Ref.  1).  In  the  previous  report,  a  system  at  a  given  temperature,  total  atmospheric  pressure, 
and  specific  humidity  is  moved  adiabatically  to  a  dry-air  pressure  of  KKX)  mb,  with  the 
resultant  temperature  of  this  new  system  being  the  potential  temperature.  Then  the  specific 
humidity  was  added  to  this  system,  and  it  was  treated  as  if  it  were  an  ideal  gas  of  the  same 
density  as  dry-air.  This  results  in  an  equivalent  dry-air  system  at  approximately  the  same 
pressure  at  a  temperature  given  by  the  virtual  potential  temperature.  This  procedure  has 
several  conceptual  and  physical  problems,  particularly  with  respect  to  the  total  atmospheric 
pressure  and  the  dry-air  pressure,  that  have  led  to  its  replacement  by  the  following  conceptual 
process.  For  a  system  at  a  given  temperature,  specific  humidity,  and  total  atmospheric 
pressure,  the  virtual  temperature  is  calculated  for  an  equivalent  dry-air  system  at  the  same 
total  atmospheric  pressure.  This  equivalent  dry-air  system  is  then  adiabatically  moved  to  a 
new  system  at  a  uniform  pressure  of  1000  mb  at  a  new  temperature  defined  as  the  virtual 
potential  temperature. 


2.2.1  Virtual  Temperature  and  Specific  Humidih 


The  virtual  temperature,  T^,  is  the  temperature  at  which  a  volume  of  dry  air  has  the  same 
density,  p,  as  moist  air  with  specific  humidity,  q,  at  the  same  total  atmospheric  pressure  (P). 
Equation  6  shows  the  relationship  between  the  virtual  temperature  and  the  temperature  of 
moist  air  (T),  and  how  this  is  applied  to  the  ideal  gas  law  (Eq.  6c): 
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[1  +  (1  -  pVp]  T 

(6a) 

MJM^  =  0.622 

(6b) 

(6c) 

In  these  equations,  is  the  molecular  mass  of  water,  is  the  molecular  mass  of  dry  air, 
is  the  gas  constant  for  moist  air,  (=  R/MJ  is  the  gas  constant  for  dry  air,  and  R  is  the 
universal  gas  constant.  In  turn,  the  specific  humidity  can  be  expressed  in  terms  of  the  total 
atmospheric  pressure  (P=P3+P,^),  the  dry  air  pressure  (P  J,  and  the  water  vapour  pressure  (P,,,) 
by 


q  =  pPy[P-(l-p)PJ 

(7a) 

(7b) 

=  (D/fl+ffl)  ;co=pPyP^  , 

(7c) 

where  ro  is  called  the  mixing  ratio.  The  new  PIRAM  model  also  follows  these  definitions. 

2.2.2  Virtual  Potential  Temperature 

The  virtual  potential  temperature  of  dry  air  (0^)  is  defined  as  the  temperature  that  a 
volume  of  dry  air  at  virtual  temperature,  T^,  assumes  when  brought  adiabatically  from  its 
existing  pressure,  P,  to  a  standard  reference  pressure,  P^^f.  The  resultant  virtual  potential 
temperature  is  related  to  the  virtual  temperature  by 

e,  "  (8) 

where  y  (=  Cp/c^  =1.40  for  dry  air  at  room  temperature)  is  the  ratio  of  the  specific  heat  at 
constant  pressure  (Cp)  and  the  specific  heat  at  constant  volume  (c^)  for  dry  air.  By  default 
and  by  convention  P„f  is  taken  to  be  1000  mb  although  L(W)WKD  does  allow  it  to  be 
modified.  In  earlier  versions  of  L(W)WKD,  the  dry-air  pressure  was  used  by  default  in  Eq.  8 
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instead  of  the  total  air  pressure.  While  L(W)WKD  still  provides  this  as  an  option,  the  total 
air  pressure  is  now  the  default  option.  The  PIRAM  model  uses  the  same  formulation. 


2.2.3  Water  Vapour  Pressure 


As  the  water  vapour  pressure  is  not  generally  available  from  meteorological  measure¬ 
ments,  L(W)WKD  calculates  it  using  the  following  empirical  relationship  for  the  saturated 
water  vapour  pressure  obtained  from  a  variation  of  the  Clausius-Clapeyron  equation  (Ref.  1 ) : 

ln(Pj(7))  =  19.32  -  4223/(r  -  32)  (9) 

Thus,  Eq.  9  provides  a  relationship  between  either  the  saturated  water  vapour  pressure  (P^^ 
and  the  air  temperature  (T),  or  using  Eq.  10b,  the  water  vapour  pressure  (P,^,)  and  the 
dew-point  temperature  (T^).  Thus,  given  the  dew-point  temperature,  the  water  vapour 
pressure  is  easily  determined.  However,  as  the  relative  humidity  (Hr)  is  more  commonly 
quoted  than  the  dew-point  temperature,  use  must  be  made  of  Eq.  10a,  which  relates  the  water 
vapour  pressure  to  the  relative  humidity,  and  the  saturated  water  vapour  pressure: 


=  H^T)  P:(T) 

(10a) 

p^{T)  =  p:(t,). 

(10b) 

The  PIRAM  model  uses  a  different  formulation  of  Eq.  9  (see  Ref.  27)  for  the  saturated  water 
vapour  pressure;  however,  this  shouldn’t  produce  any  significant  differences. 


At  the  surface  of  a  sea  with  salinity,  S,  and  at  temperature  T^,  the  relative  humidity  is  not 
100%  as  was  previously  used  by  default  in  most  earlier  versions  of  both  L(W)WKD  and 
PIRAM.  Instead,  it  is  given  by 


PL(T,.S)  ^  p:(TJ(1  -  eS) 

p:(tj  p:(tj 


=  (1  -eS) 


(lla,b) 
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where  is  the  saturated  water  vapour  pressure  above  salt  water  (sw),  is  the  saturated 
water  vapour  pressure  above  fresh  water,  and  e  is  taken  to  be  5.37  x  10'^  (Ref.  13).  Thus,  for 
a  standard  salinity  of  34g  of  salt  per  kilogram  of  water,  the  relative  humidity  at  the  surface 
is  98.2  %.  This  is  now  the  default  value  for  the  relative  humidity  at  the  sea  surface  for  both 
L(W)WKD  and  PIRAM. 


2.2.4  Vertical  Atmospheric  Pressure  Profile 


The  final  required  equation  is  that  relating  the  total  air  pressure  (P  =  P^  +  P^)  to  the 
height  above  the  ground.  These  equations  are  different  from  those  in  the  earlier  versions  of 
L(W)WKD  due  to  the  change  in  conception  of  the  physical  processes  discussed  earlier. 
Assuming  that  the  forces  on  an  equivalent  slab  of  dry-air  are  balanced,  and  using  Eq.  6c  for 
the  perfect  gas  law,  gives  the  hydrostatic  equation: 


^  -  -_A_ 
P  KT. 


dz. 


(12) 


where  dP  is  the  pressure  differential  across  the  slab,  dz  is  the  width  of  the  slab,  g  is  the 
acceleration  due  to  gravity  near  the  surface,  T^  is  the  virtual  temperature  of  the  slab,  and  R^ 
is  the  gas  constant  for  dry-air.  Assuming  that  g  and  R^  do  not  vary  much  near  the  surface, 
integrating  Eq.  12  gives; 


P(z)/P{z,) 


S  r  dz' 


-g(z-zi) 


j  _  g(z-Zi) 


(13a,b,c) 


where  P(Zi)  is  the  total  atmospheric  pressure,  and  Ty(Zi)  is  the  virtual  air  temperature  at  the 
height  z,.  The  approximation,  Eq.  13b,  is  used  by  L(W)WKD,  and  is  obtained  by  assuming 
that  the  temperature  profile  is  uniform  (constant  temperature  within  the  MBL).  The  further 
simplification  of  Eq.  13c  is  obtained  by  using  the  first  term  of  the  Taylor  series  expansion 
of  the  exponential  in  Eq.  13b.  This  final  form  is  used  in  PIRAM;  however,  it  should  not 
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produce  any  significant  differences  between  the  models.  The  air  density  profile  can  be 
obtained  from  Eq.  13b  and  Eq.  6c,  such  that 


-g(z-Zl) 

p(z)  =  -^e 


(14) 


Again  PIRAM  expands  the  exponential  and  only  uses  the  first  term  in  z. 


2.3  Solving  the  Basic  Scaling  Equations 

To  solve  Eq.  4  requires  the  determination  of  ten  unknown  parameters.  Six  of  these 
parameters  can  be  determined  using  boundary  conditions,  while  the  other  four  need  to  be 
determined  using  other  techniques.  The  current  differences  between  the  previous  L(W)WKD 
versions  and  the  PIRAM  model  are  highlighted  in  the  following  sections. 


2.3.1  Boundary  Conditions 


At  sea  level,  the  following  three  boundary  conditions  allow  the  following  definitions  to 
be  made  for  the  three  roughness  lengths.  First,  let  be  the  height  at  which  the  wind  speed 
is  zero  (u(Zou)  =  0),  Zpq  the  height  at  which  the  water  vapour  pressure  is  equal  to  the  sea 
surface  water  vapour  pressure,  and  Zqq  the  height  at  which  the  virtual  potential  temperature 
of  the  air  is  equal  to  that  given  by  air  at  the  temperature  of  the  sea  (T^),  and  the  water  vapour 
pressure  is  given  by  the  sea  surface  water  vapour  pressure.  Applying  these  definitions  to 
Eq.  4,  and  using  Eqs.  6,  7,  and  8  gives: 

“o  "  ®  (15a) 

q,=q{z,^)  =  p//,(r;pJ(r;/[P,(0)+M//,(r;pJ(r;]  (isb) 

+(1  -P) (15c) 


where  Pa(z=0)  is  the  dry-air  atmospheric  pressure  in  millibars,  P(z=0)  is  the  total  atmospheric 
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pressure,  HrCTJ  is  the  relative  humidity,  and  P^'CTj)  is  the  saturated  water  vapour  pressure 
at  the  water  level  for  water  at  temperature,  T^. 

Three  more  boundary  conditions  are  given  using  measurements  of  the  wind  speed  u(Zo) 
at  the  height  Zq,  and  measurements  of  the  air  temperature,  T(Zi),  total  atmospheric  pressure, 
P(Zi),  and  relative  humidity,  H,(Zi)  at  the  height,  z,.  Applying  these  conditions  to  Eq.  4,  and 
using  Eqs.  6, 7  and  8  gives: 

Au(Zo)=  u{Zq)  =(M,/K)[ln(Z(/ZoJ+T„(Zo;L)-T„(Zo„;L)]  (16a) 

Ae,(Zi)  =e,(Zi)  -0^0  =  (0vA)[ln(Zi/Zo0,)  "^(^00/^)]  (16b) 

Aq{z^)  =  q{z^) -qQ  =(^yK)[ln(Zj/ZoP  (I6c) 

where 

q{z,)  =MH,a,)Pj(7’(2i))/[P,(z,)  (17a) 

ev(2i)  =  7’(^i)[l  +(1  -p)^(z,)/p][l{)00/P(z,)r->>^  (17b) 

q(Zi )  is  the  specific  humidity,  0^(zi)  is  the  virtual  potential  temperature,  and  P^®  (T(Zi))  is 
the  saturated  water  vapour  pressure  at  height  Zj.  Furthermore,  A0v(Zi)  is  known  as  the  virtual 
potential  air-sea  temperature  difference  (VPASTD),  and  Aq(Zi)  is  known  as  the  air-sea 
specific  humidity  difference  (ASSH). 

The  two  sets  of  boundary  conditions  provide  six  equations  for  the  determination  of  the 
10  unknowns  and  are  used  by  both  L(W)WKD  and  PIRAM.  The  following  two  sections 
provide  the  remaining  four  equations. 

2.3.2  Roughness  Lengths 


The  following  empirical  relationships  have  been  developed  by  Businger  et.  al.  (Ref.  14) 
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and  Smith  (Ref.  15).  They  relate  to  Zqu,  and  2^  to  by; 

ln(zgg/10)  =  -K^/lCj^lnClO/zg^)]  (18a) 

,  ln(Zo/10)  =  -K2/[C^^ln(10/Zo„)]  (18b) 


where  Cjn  is  the  10  m  neutral  heat  flux  coefficient  whose  value  has  been  established  to  be 
1.0  X  10'^  for  wind  speeds  ranging  from  6  to  22  m/s,  and  Cen  is  the  10  m  moisture  flux 
coefficient  defined  by  Smith  (Ref.  15)  and  evaluated  to  be  1.2  x  10'^.  To  solve  Eq.  18,  use 
is  made  of  the  following  wind  dependent  relation  (Ref.  1,9,10  &  15); 


z 


OU 


{ajg)ul  + 
(aJg)K  + 


V^-CK 


e 


-Ck 


(19a) 

(19b) 


where  k  is  von  Karmen’s  constant,  C  is  a  constant  equal  to  -5.5,  g  is  the  acceleration  due  to 
gravity,  (=  0.01 1)  is  Chamock’s  constant,  and  v^  is  the  dynamic  viscosity  for  diy  air.  As 
with  the  previous  versions  of  L(W)WKD,  these  values  have  been  chosen  according  to 
Walmsley  (Ref.  1 1).  The  difference  with  the  earlier  versions  is  that  the  dynamic  viscosity 
is  now  taken  to  be  a  function  of  the  virtual  potential  temperature  (K),  where  is  9.267  x  10'® 
and  is  1.346  x  10'^  (see  Eq.  19b).  These  values  are  obtained  from  a  linear  fit  of  the  data 
given  in  Table  A3  of  Ref.  16  for  temperatures  from  0  to  40°C  at  a  pressure  of  1000  mb.  Both 
L(W)WKD  and  PIRAM  now  use  this  functional  form.  However,  L(W)WKD  uses  the 
average  virtual  potential  temperature,  given  by; 

+  0v(^,)]/2  =  0,  ,  (20) 

where  is  the  virtual  potential  temperature  at  the  sea  surface,  and  0v(Zi)  is  the  virtual 
potential  temperature  at  height  Zj,  in  place  of  the  virtual  potential  temperature  given  in  Eq. 
19b.  The  exact  formulation  used  by  PIRAM  is  not  known;  however,  it  is  doubtful  that  any 
significant  differences  would  result  in  a  slightly  different  interpretation. 
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2.3.3  Monin-Obukhov  Length 

With  the  addition  of  Eqs.  18  and  19  there  are  now  9  equations  for  10  unknowns  with  the 
remaining  unknown  being  the  Monin-Obukhov  length  (L).  Equation  2  provides  the  required 
10th  equation  provided  a  definition  is  given  for  the  average  virtual  temperature,  As  in 
previous  versions  of  L(W)WKD,  this  average  is  taken  between  the  respective  temperature 
at  the  sea  surface,  and  that  measured  at  the  height  z^; 

K.  =  IKs  "  ^v(^i)]/2.  (21) 

Here,  again  as  for  Eq.  20,  L(W)WKD  and  PIRAM  may  have  slightly  different  interpretations; 
however,  again,  no  significant  differences  are  expected. 

2.4  Iterative  Algorithm 

The  algorithm  described  below  is  the  one  used  by  the  new  L(W)WKD  program  to 
provide  solutions  for  any  set  of  atmospheric  parameters  of  interest  and  has  been  slightly 
modified  from  the  earlier  versions  (before  Ver.  7.0).  The  algorithm  is  enumerated  as  a  series 
of  steps,  with  each  step  accompanied  by  a  descriptive  comment.  For  simplicity  the  algorithm 
is  divided  into  three  parts.  The  first  part  describes  the  steps  required  to  determine  the 
measured  virtual  potential  air-sea  temperature  difference  (VPASTD  =  AG^)  and  the  air-sea 
specific  humidity  difference  (ASSHD  =  Aq),  while  the  second  part  describes  the  steps  used 
to  initialize  the  variables  used  during  the  iteration  algorithm  described  in  the  third  section. 
The  iterative  algorithms  used  by  L(W)WKD  and  PIRAM  still  remain  different  as  they  have 
not  been  part  of  any  discussions  to  bring  them  into  agreement. 

2.4.1  Part  1  -  Measurements 


1)  Calculate  the  water  vapour  pressure,  using  Eqs.  9  and  10a  with  the 

measured  relative  humidity  and  air  temperature  at  the  height  Zj. 
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2)  Calculate  the  specific  humidity,  qfzj),  using  Eq.  7a  with  the  measured  air 
pressure  at  the  height  Zj,  and  the  water  vapour  pressure  from  Step  1. 

3)  Calculate  the  virtual  temperature,  T^Czi)  using  Eq.  6a,  the  measured  air 
temperature  at  the  height  Zj,  and  the  specific  humidity  from  Step  2. 

4)  Calculate  the  virtual  potential  temperature,  0v(zi)  using  Eq.  8,  the  measured  air 
pressure  at  the  height  Zj,  P(Zi),  and  the  virtual  temperature  from  Step  3.  If  the  dry- 
air  option  has  been  chosen  then  the  air  pressure  is  replaced  by  P(Zi)  -  Pw(Zi),  where 
P,,,(Zi)  is  obtained  from  Step  1. 

5)  Calculate  the  total  atmospheric  pressure  at  sea  level,  P(z=0),  using  Eq.  13b, 
the  total  air  pressure  at  height  z,,  and  the  virtual  temperature  from  Step  3. 

6)  Calculate  the  water  vapour  pressure  at  sea  level,  P,^,(z~0)  =  P\„(z!=0),  using 
Eqs.  9  and  1  la  with  the  sea  surface  relative  humidity  and  sea  surface  temperature. 

7)  Calculate  the  specific  humidity  at  sea  level,  q(z«0),  using  Eq.  7a  with  the  air 
pressure  from  Step  5  and  the  water  vapour  pressure  from  Step  6. 

8)  Calculate  the  virtual  temperature  at  sea  level,  Ty(z»0)  using  Eq.  6a,  the 
measured  sea  temperature  and  the  specific  humidity  from  Step  7. 

9)  Calculate  the  virtual  potential  temperature,  9v(z=0)  using  Eq.  8,  the  air 
pressure  from  Step  5,  and  the  virtual  temperature  from  Step  8.  If  the  dry-air  option 
has  been  chosen  then  the  air  pressure  is  replaced  by  P(z«0)  -  P,y(z«0),  where 
P,,,(z~0)  is  the  result  from  Step  6. 

10)  Calculate  the  average  virtual  temperature  in  the  layer,  T^,  using  Eq.  21  and 
the  results  from  Steps  3  and  8. 


11)  Calculate  the  average  virtual  potential  temperature  in  the  layer,  using 
Eq.  20  and  the  results  from  Steps  4  and  9. 


2.4.2  Part  2  -  Initialization 


1)  Set  the  number  of  iterations  parameter,  N,  equal  to  zero.  This  parameter  keeps 
track  of  the  number  of  iterations  required  to  obtain  a  stable  solution. 
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2)  Set  the  iteration  precision  parameter,  e,  equal  to  0.00 1 .  This  is  one  of  the  two 
arbitrarily  set  parameters  that  can  stop  the  iteration  procedure.  Its  significance  is 
discussed  further  in  Step  10  of  Part  3. 

3)  Set  the  maximum  number  of  iterations  parameter,  MAXN,  equal  to  1 00.  This 
is  the  second  arbitrary  set  parameter  that  can  stop  the  iteration  procedure. 

4)  Set  the  following  parameter,  S(0),  equal  to  a  large  number.  It  is  currently  set 
arbitrarily  to  1  x  10’”.  The  significance  of  this  parameter  is  discussed  in  Step  10 
of  Part  3. 

5)  Calculate  the  air-sea  specific  humidity  difference  (ASSHD  =  q(Zi)  -  qo)  using 
the  results  from  Steps  2  and  7  of  Part  1. 

6)  Calculate  the  virtual  potential  air-sea  temperature  difference  (VPASTD  = 
0,,(z,)  -  0^)  using  the  results  from  Steps  4  and  9  of  Part  1. 

7)  Set  the  wind  speed  measured  at  the  height  Zq  to  u(Zo). 

8)  Set  the  initial  wind  roughness  height,  Zou(O)  equal  to  1.5  x  lO"^  m.  This  is  a 
generally  accepted  approximate  value  (Ref.  1). 

9)  Set  the  initial  virtual  potential  temperature  roughness  height,  Zoev(O),  using 
Eq.  18a. 

10)  Set  the  initial  specific  humidity  roughness  height,  Zoq(O),  using  Eq.  18b. 

11)  Set  the  initial  wind  scaling  constant,  u*(0),  using  only  the  first  term  in  Eq.  19 
and  the  value  of  Zqu  from  Step  8: 

M,  =  [zJ\.\2xm-Y^  =  0.3660  mis.  (22) 

12)  Set  the  initial  virtual  potential  temperamre  scaling  constant,  0^.(0),  to 
(k  X  VPASTD)/z,  where  k  is  von  Karmen’s  constant  and  VPASTD  is  obtained 
from  Step  6. 

13)  Set  the  initial  specific  humidity  scaling  constant,  q.(0),  to  (k  x  ASSHD)/Zi 
where  k  is  von  Karmen’s  constant  and  ASSHD  is  obtained  from  Step  7. 

14)  Set  the  initial  Monin-Obukhov  length,  L(0),  using  Eq.  2  and  the  results  from 
Steps  11  and  12,  and  Step  10  from  Part  1. 
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2.4.3  Part  3  -  Iteration 


1)  Using  the  wind  speed  roughness  height,  Zou(O),  the  Monin-Obukhov  length, 
L(0),  and  u(Zo),  solve  Eq.  16a  for  a  new  value  of  the  scaling  constant,  u.(l). 

2)  Using  the  virtual  potential  temperature  roughness  height,  (0) ,  the  Monin- 
Obukhov  length,  L(0),  and  the  VPASTD,  solve  Eq.  16b  for  a  new  value  of  the 
scaling  constant,  0v.(l). 

3)  Using  the  specific  humidity  roughness  height,  Zq(0),  the  Monin-Obukhov 
length,  L(0),  and  the  ASSHD,  solve  Eq.  16c  for  a  new  value  of  the  scaling  constant, 

q*(l). 


4)  Calculate  a  new  value  for  the  wind  speed  roughness  height,  Zou(l),  using 
Eq.  19b,  u*(l)  from  Step  1,  and  from  Step  11  of  Part  1. 

5)  Calculate  a  new  value  for  the  virtual  potential  temperature  roughness  height. 
Zoo  (1),  using  Eq.  18a  and  Zou(l)  from  Step  4. 

6)  Calculate  a  new  value  for  the  specific  humidity  roughness  height,  ZgqC  1 ),  using 
Eq.  18b  and  Zou(l)  from  Step  4. 

7)  Calculate  a  new  value  for  the  Monin-Obukhov  length,  L(l),  using  Eq.  2, 
from  Step  10  of  Part  1,  and  the  new  scaling  constants  u*(l)  and  0v»(l)  calculated 
from  Steps  1  and  2. 

8)  Calculate  the  absolute  relative  differences  for  the  three  scaling  constants  and 
the  three  roughness  heights.  The  absolute  relative  difference  is  defined  to  be: 


D, 


j[A(l)  -  A(0)]/A(0)1  ;  A(0)^0 
‘|[A(1)  -A(0)]/A(l)i  ;  A(0)=0,A(1)^0 
0  ;A(0)=A(1)=0 


(23) 


where  A  represents  any  of  the  six  parameters. 

9)  Calculate  the  sum,  S(l),  of  the  six  absolute  relative  differences  calculated  in 
Step  8.  S  is  the  measure  used  to  describe  the  goodness  of  the  solution. 
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10)  The  final  step  in  the  iteration  loop  is  the  following  IF  block  which  decides 
whether  to  continue  or  stop  the  iterative  procedure.  The  procedure  is  stopped  and 
the  solution  is  valid  if  S(l)  is  less  than  s  and  N  is  less  than  MAXN.  The 
significance  of  S(l)  being  less  than  s  is  that  the  sum  of  the  relative  differences  is 
less  than  100*8  %  (i.e.  the  maximum  possible  value  that  any  of  the  individual 
relative  differences  can  have  is  100*8  %).  Thus,  a  value  of  0.001  translates  into  a 
maximum  difference  of  0.1%  for  any  of  the  parameters.  The  IF  block  used  by 
L(W)WICD  is  constructed  as  follows: 

if(S(1)<8)  then 

-  solution  found;  stop  the  iteration  procedure  and  return  the  results 
ELSE  IF  (N  >  MAXN)  THEN 

-  no  solution  found  within  MAXN  iterations;  stop  the  iteration  procedure 
and  return  an  error  flag 

ELSE  IF  ((N  ^  5)  OR  (S(l)  <  S(0))  THEN 

-  solution  not  yet  found;  continue  the  iteration  procedure  as  either  less 
than  6  iterations  have  been  performed  (N  s  5)  or  the  goodness  of  the 
solution  is  improving  (S(l)  <  S(0)) 

-  increment  N  to  N+1 

-  set  Zou(O)  =  zjl);  ZosvfO)  =  Zoev(l);  Zoq(O)  =  Zoq(l) 

-  set  u.(0)  =  u.(l);  0,.(O)  =  0,.(1);  q.(0)  =  q.(l) 

-  set  L(0)  =  L(l) 

-  GOTO  Step  1 

ELSE 

-  no  solution  found  because  N  >  5  and  S(l)  ^  S(0) 

IF  (mode  is  INTERACTIVE*)  THEN 

-  stop  the  iteration  procedure  and  return  an  error  flag 
ELSE  (mode  is  BATCH*) 

IF  (wind  speed  <  25  m/s)  THEN 

-  increase  the  wind  speed  by  0.5  m/s  and  return  to  Step  1 
ELSE 

-  stop  the  iteration  procedure,  return  an  error  flag,  and 
output  uniform  profiles 


*  The  INTERACTIVE  and  BATCH  modes  of  operation  are  discussed  in  Section  3. 
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2.5  Waw  Model 

The  Wavy  WKD  model  developed  at  DREV  by  Beaulieu  (Ref.  10)  was  an  attempt  to 
mimic  the  effect  that  water  waves  would  have  on  the  vertical  profiles  generated  by  the  WKD 
model  (flat  sea).  The  L(W)WKD  model  also  incorporates  this  effect  on  the  profiles  it 
generates.  This  feature  is  not  present  in  PIRAM  and  will  only  be  discussed  briefly  in  this 
report  as  Ref.  1  discusses  it  in  detail. 


If  one  assumes  that  the  action  of  the  waves  is  to  turn  a  flat  surface  into  a  sinusoidal 
surface  with  a  wave  height  (trough  to  peak)  of  and  a  wavelength  of  then  the  height  of 
any  part  of  the  wave  (H)  above  the  mean  sea  level  is  given  by 

H(x)  =  iHJ2)cos(2TzxrkJ,  (24) 


where  x  is  between  0  and  Thus,  assuming  that  the  instantaneous  vertical  profile,  0„,  for 
parameter  a,  above  every  point  of  the  wave  is  that  given  by  L(W)WKD  for  a  flat  sea,  the 
average  effect  of  this  wave  action  can  be  calculated  by  averaging  the  vertical  profile  over  one 
wavelength.  This  gives: 


K(z)  = 


Z*HJ2 


iHHJ  f  0^(z  )dz '  ;z^z^+HJ2 


Z-HJ2 


z^HJ2 


(25) 


( HHJ  J  0„(z ')dz'  ■,z^<z< z^^  +  HJ2 


^inax 


where  O/  is  the  averaged  vertical  a-profile,  z  is  the  height  above  the  mean  sea  level  or  flat 
sea,  z'  is  the  instantaneous  height  above  a  point  on  the  wave,  and  z^^  (=max(Zou,  z^^,  z^) 
is  the  maximum  of  the  wind  speed,  virtual  potential  temperature,  and  specific  humidity 
roughness  heights.  The  first  function  in  Eq.  25  is  used  for  heights  above  the  peak  of  the 
wave,  while  the  second  function  is  used  for  heights  between  the  mean  sea  level  and  the 
wave’s  peak.  The  function  is  not  calculated  for  heights  below  the  mean  sea  level. 
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Equation  25  can  easily  be  calculated  using  the  rectangular  integration  rule  such  that 

n-l 

E HJ1\  i  =0,1, ..n-l 


Kiz)  = 


(=0 

n~l 


(26) 


E ^„(z/)/(« -^z^<z<z^+  HJ2-  i  =0, 1  ,..n - 1 


i=m 


z/  =  z~iHJl)cos{lniln) 


(27) 


and 

m  =  (n/3i)cos-*[2(z-z^J///J.  (28) 

The  above  procedure  is  only  valid  for  profiles  that  are  linear  combinations  of  the  four 
fundamental  vertical  profiles  given  by  Eqs.  4a,  4b,  4c  and  13b  for  the  wind,  virtual  potential 
temperature,  specific  humidity,  and  total  atmospheric  pressure  respectively.  Any  nonlinear 
combinations  of  these  profiles  must  first  calculate  the  wavy  profile  for  each  of  these  four 
profiles  before  combining  their  results.  This  is  an  extremely  important  point  (see  Ref  1),  as 
many  of  the  desired  quantities,  involve  nonlinear  combinations  of  these  profiles. 

2.6  Refractive  Index 


Discussion  of  the  functional  form  for  the  index  of  refraction  or  refractivity  will  be 
relatively  brief  as  it  is  fully  discussed  in  Ref  1.  While  the  functional  form  used  by 
L(W)WKD  and  PIRAM  is  the  same,  the  same  is  not  true  for  the  various  coefficients  and  this 
has  not  been  part  of  any  common  agreement. 

Following  Ref  1,  the  refractivity  (Z  =  (n-l)  x  10®,  where  n  is  the  refractive  index),  of 
the  atmosphere  at  a  wavelength,  X,  can  be  approximated  by 
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Z(k)x  10-^  =  \D,(KT)^  + 


(100  Pa/mb) 


P  P 

A(k,T)-^  +  B(k,T)^\  X  10'^ 


(29a) 

(29b) 


where  k  is  Boltzmann’s  constant,  Z  is  the  refractivity,  is  the  dry-air  pressure,  is  the 
water  vapour  pressure,  T  is  the  temperature,  is  the  dry-air  molecular  refractivity,  is 
the  water  vapour  molecular  refractivity,  A  is  the  dry-air  coefficient,  and  B  is  the  water  vapour 
coefficient.  In  all  future  discussions  it  will  be  assumed  that  all  pressures  are  in  millibars,  all 
temperatures  are  in  Kelvin,  all  densities  are  in  grams  per  cubic  metre,  and  all  wavelengths 
are  in  micrometers.  This  leads  to  various  transformation  constants  in  many  of  the  equations, 
including  the  factor  of  100  Pascal/1  millibar  given  in  Eq.  29a. 


2.6.1  Dry- Air  Refractivity 


Using  measured  refractivities  for  standard  air  (defined  as  containing  0.03%  COj  per 
volume  at  101.3  mb  and  15  °C),  many  authors  have  determined  the  A  coefficient  of  Eq.  29. 
As  discussed  in  Ref.  1,  L(W)WKD  uses  the  following  Sellmeier  function  (Ref.  17) 
developed  by  Edlen  in  1953  (Ref.  18)  for  A: 


A(K)  = 


64.328 


29498.10  255.40 

+  + 

288.2 

146  - 1/^2  41  -  1/k^ 

1013 

(30a) 


=  a(k) 


288.2 

1013 


(30b) 


where 


a(X,) 


64.328. ^94^ 

146-1/^2 


(31) 


k  is  the  wavelength  in  vacuum  in  pm  and  k  is  Boltzmann’s  constant.  It  has  simple  poles  at 
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0.083  |jm  (=[1/146]*^)  and  0.156  |am  (=[1/41]*'^)  which  simulate  the  effects  of  resonances 
due  to  bound  states,  states  near  the  ionization  limit,  and  dissociative  states.  Edlen  has  shown 
that  this  fit  is  good  to  within  0.04  units  (~  .05%)  down  to  0.2  ^m. 


2.6.2  Water  Vapour  Refractivitv 


For  water  vapour,  L(W)WKD  uses  the  results  of  Hill  and  Lawrence  (Ref.  7),  which 
separates  the  molecular  refractivity  into  the  following  sums  over  the  IR  and  UV  resonances 
(?^); 

IR  UV 

(32) 

=Dl\l,T)  +  D^iX,T)  , 


where  Sellmeier's  dispersion  function  (Ref.  17)  is  used,  and  the  Sj’s  are  the  line  strengths. 
Using  dispersion  data  from  Erickson  (Ref.  19)  for  eight  visible  wavelengths  between  0.644 
and  0.361  pm,  they  give  the  UV  molecular  refractivity  to  be 


-  37.47  k 

124.40 


(33a) 


where 


= 


31A1 

124.40 -l/X^  ’ 


(33b) 


and  R^  (=R/M,j,=0.4619)  is  the  gas  constant  for  water  vapour.  There  is  a  simple  pole  at 
0.090  pm  which  simulates  the  resonances  from  bound  states,  states  near  the  ionization  limit, 
and  dissociation  states.  As  for  dry-air,  this  implies  that  Eq.  33  should  be  good  down  to 
0.2  pm.  The  claimed  accuracy  is  less  than  0.5%.  A  plot  of  b^(X.)  is  given  in  Fig.  1. 
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UV  &  IR  b  Coefficient  vs  Wavelength 


□  UV  A  IR;-40C  A  IR;-hWC 
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FIGURE  1  -  Plot  of  b  coefficients  (x  100/R^  =  216.7) 


To  determine  the  refractivity  due  to  the  IR  resonances,  Hill  and  Lawrence  (Ref.  7) 
calculated  D„®(X,,T)  by  performing  the  summation  over  more  than  40,000  water-vapour 
resonances  extending  from  radio  frequencies  to  the  visible  for  nine  wavelengths  and 
temperatures  from  -40  to  40  ”C  in  10  °C  steps.  The  wavelengths  were  18.83, 17.82, 15.43, 
13.1 1,  1 1.60,  10.37,  9.341,  8.562  and  7.837  pm.  After  fitting  the  calculations  and  scaling 
the  results  to  agree  with  measurements  near  10.6  pm,  the  following  IR  refractivity  function 
was  found. 


dI\X,T)  = 


0.957 -0.928  6°  \jc-l) 


1.036°  -  19.8x2  +  8.1jc^  -  1.7x^ 


JL=b^R(k,T)  — 

K  K 


(34a) 


where 
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b'%,T)  = 


0.957- 0.928  5°  \;c-l) 
1.036°  -  19.8>:2  +  8.1a:^  - 1.7  ;c^ 


(34b) 


6  =  T/273. 15,  and  x  =  10/X.  The  plot  of  b®(A,,T),  in  Fig.  1,  shows  the  effect  of  changing  the 
temperature.  As  can  be  seen,  the  effect  of  temperature  is  not  very  important,  so  that  using 
8=1  does  not  significantly  affect  the  output  of  Eq.  34.  This  graph  is  different  from  that 
given  in  the  previous  report  (Ref.  1)  due  to  an  error  of  sign  that  was  introduced  when  creating 
the  graph.  Fortunately,  this  error  does  not  produce  a  significant  effect  on  the  values  of  the 
A  and  B  coefficients  (see  Table  1)  in  the  regions  of  interest.  Adding  the  UV  and  IR 
contributions  together,  we  get  the  following  result  for  the  refractivity  of  water  vapour; 

Uy  jf,  IXIO^P  P 

ZJ}.,T)  =  [D^(k)^D%,T)-\~- .  _  ^  (35a) 

K1  X 

where 

R(?.,7)xl0-«  =  =  [b^(X)+b’\X,T)]/R^  =  b(X)/R^  ,  (35b) 

and 

b(k)  =  [b^(k)+b^‘^(k,T)]  .  (35b) 

Equation  35  is  applicable  for  wavelengths  between  0.2  and  19  pm  when  they  are  sufficiently 
far  from  any  significant  resonances.  In  particular,  this  means  that  Eq.  35  should  work  very 
well  in  the  visible,  the  3-5  pm  and  the  8-12  pm  windows. 

2.6.4  The  A  &  B  Coefficients 


Figure  2  shows  a  plot  of  the  A  and  B  coefficients  defined  in  the  preceding  sections  for 
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A,  B  Coefficient  vs  Wavelength 


[-g-A  (Edien  53)  -^-B  (Hill  (86) 


FIGURE  2  -  Plot  of  A  and  B  coefficients  versus  ivavelength 
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0  °C;  however,  the  L(W)WKD  program  does  not  compute  these  coefficients  for  any  desired 

wavelength.  Instead,  it  uses  precomputed  averages  over  both  temperature  and  wavelength 

for  four  bands  between  0.4  and  12  pm,  and  for  the  microwave/radar  band.  Table  I  shows  the 

values  for  A  and  B  along  with  a  measure  of  their  uncertainty.  These  values  are  essentially 

those  given  by  Ref.  1  except  for  some  minor  changes  in  the  B  coefficient.  For  ease  of 

comparison,  its  values  from  Ref.  1  are  given  in  the  last  column.  The  PIRAM  model  is  very 

similar  except  that  it  does  not  use  band  averages  for  the  visible  and  IR.  Instead,  it  calculates 

the  A  coefficient  for  a  specified  wavelength,  and  the  B  coefficient  from  the  A  coefficient 

using  B  =  A  -  12.79.  As  can  be  seen  from  Table  1,  this  does  not  cause  any  significant 

differences  except  in  the  far-infrared.  The  values  used  in  the  RF  band  are  taken  from 
Ref.  20. 
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TABLE  I 


A  &  B  Wavelength  Coefficients 


Wavelength  Band 

XCgm) 

B  (K/mb)  (Ref.  1) 

Visible 

P 

1 

d 

78.6  ±  0.6 

66.6  ±  0.8 

67.0  ±  0.8 

Near-IR 

1.0 -3.0 

77.6  ±0.1 

65.5  ±0.1 

65.5  ±0.1 

Mid-m 

3.0 -5.0 

77.52  ±  0.01 

65.33  ±0.01 

65.0  ±0.2 

Far-m 

8.0  - 12.0 

77.491  ±0.001 

48.5  ±  5.4 

40.0  ±7.7 

RF  &  MicroWave 

<  50  GHz 

77.6 

72  +  3.75  X  lO’/T 

3.734  X  10^/T 

2-7  Structure  Parameters 


The  L(W)WKD  model  also  calculates  profiles  for  the  refractive  index  structure 
parameter,  C„2,  the  refractivity  index  structure  parameter,  Cz2,  and  other  structure  parameters. 
Differences  with  PIRAM  continue  to  exist  in  this  area  as  it  was  not  part  of  our  common 
agreement.  As  such,  this  section  is  relatively  brief  as  the  reader  is  referred  to  Ref.  1  for  more 
detailed  information. 

As  shown  in  Ref.  1,  the  refractive  index  structure  parameter  at  an  elevation,  z,  is  given 
by 

C„2(z)  =  h~'^^<[n{z)-n{z+h)Y>^^ 

(36) 

=  1x10-12  h-^^<[Z{z)-Z{zm\,,  -  lxlO-^^C^2(z) 

where  h  is  some  incremental  height,  Z  is  the  refractivity,  and  Cz2  is  the  refractivity  structure 
parameter.  Using  Eq.  29b,  the  refractivity  structure  parameter  can  be  expressed  as 

h^^C^2(z) -<[{K^(zXT)iT(z)-T(z^h))  ^K^a,T)(p^(z)-pXz^h))]\^^  (37) 


where 
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K^izXT) 


A(k)P,iz) 

T\z) 


Y.JzX,T)  =  BXT) 


R. 


100  Paimb 


(38a, b) 


Carrying  out  the  square  in  Eq.  37,  gives 

h^^C^2{z)-<  kJ  [T{z)-T{z*h)f  +  kI  [p,(z)-p,(z+/i)]" 

mz)-T{z^h)][p^{z)-pX^m> 


(39) 


and  finally 

qafe)  =  Cr*)  ♦  kI  C^i(z)  *  .  (40) 

where  C7.2  is  the  temperature  structure  parameter,  C„  2  is  the  water  vapour  density  structure 

i  Pyv 

parameter,  and  Cj.  is  the  temperature-water  vapour  density  structure  parameter. 

•*r  W 

Following  the  parametrization  of  Monin-Obukhov,  the  virtual  potential  temperature 
structure  parameter,  the  specific  humidity  structure  parameter,  and  the  virtual  potential 
temperature-specific  humidity  structure  parameter  can  be  expressed  in  terms  of  the  scaling 
parameters,  0^*  and  q.,  such  that 


''v 

(41a) 

Cg2(z)  =  ql  Z  f^izJL) 

(41b) 

"  '■Ov?  <1*  z-^^UzJL) 

(41c) 

where  is  the  virtual  potential  temperature-specific  humidity  correlation  coefficient,  and 
the  fj’s  are  experimentally  determined  functions.  As  the  water  vapour  density  is  related  to  the 
specific  humidity  by  p,^  =  pq,  where  p  is  the  total  air  density,  we  can  express  the  water  vapour 
density  structure  parameter,  and  the  virtual  potential  temperature-water  vapour  density 
structure  parameter  by 
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Cp2(z)  =  p\z)  C^2iz)  (42a) 

C^^iz)  =  p(z)  C,^(z),  (44b) 

where  we  have  assumed  that  the  total  density  changes  very  little  between  z  and  z+h. 
Similarly,  using  Eqs.  6a  and  8,  and  assuming  that  the  specific  humidity  and  the  pressure 
change  very  little  between  z  and  z+h,  the  temperature  structure  parameter  and  the  virtual 
temperature  structure  parameter  are  given  by 

Cj.2iz)  =  C^2(z)/[1  +  (1  -  p)9(z)/p]2 

V 

CUz)  =  (P(z)/1000)2(^-^)'^  CM). 

'  V 

These  results  also  imply  the  following  two  relations: 

C,,(z)  =  Cj-^bW  +  (1  -  M)9fe)/M] 

C,/z)  =  (nz)/iooo)»-'»c,^fe) . 

Subsequent  substitution  of  the  results  from  Eqs.  41  thru  44  into  Eq.  40  allows  one  to  easily 
calculate  either  the  refractive  index  structure  parameter,  or  the  refractivity  structure 
parameter  once  we  have  some  knowledge  about  the  correlation  coefficient,  r„  ,  and  the 
structure  parameter  functions,  f/s,  in  Eq.  41. 

2.7.1  Structure  Parameter  Functions 

Several  expressions  for  the  structure  parameter  function,  fj,  have  been  proposed.  Edson 
et  al.  (Ref.  2 1 )  proposed  a  modified  version  of  that  given  by  Wyngaard  et  al.  (Ref.  22),  while 
F.F.O.*  (see  Ref.  23)  has  proposed  an  independent  function.  These  three  functions  are  given 


(43a) 

(43b) 

(44a) 

(44b) 


‘Forschuninstitut  Fiir  Optik  (Germany) 
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below: 


(Wyngaard) 


f.izJL)  = 


4.9  [1  -7(z/L)]-2/3  .  ^^<0 
4.9[l+2.75(z/L)]  ;  z/L^O 


(45) 


(Edson) 


5.8[l-7(z/I)]'^^  ;  z/L<0 
5.8  [1  +2.4(z/L)2^^]  ;  z/L^O 


(46) 


(F.F.O.)  f,(z/L)  = 


f6.3[l-7(z/L)+75(z/L)2]-i'3  ;  z/L<0 
[6.3  [  1  +  7(z/L)  +  20(z/L)2] ;  z/L  ^  0  . 


(47) 


L(W)WKD  uses  that  of  Edson  et  al.;  however,  the  other  two  functional  forms  may  be  equally 
valid.  The  PIRAM  model  uses  the  F.F.O.  functional  form. 


As  it  is  difficult  to  determine  the  function  fj,  most  authors  have  worked  under  the 
assumption  that  it  is  proportional  to  fj.  This  hypothesis  implies  that 


Cq2 


(48) 


where  the  structure  parameter  ratio  A2i(=f2/fi),  is  the  proportionality  constant  to  be 
determined,  and  which  must  always  be  positive.  Similarly,  using  the  following  definition 


(49) 


it  follows  that 


/j(z/L)  =  /2(z/L)M2,  -  f,(z/L)/Aii\ 


(50) 
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The  correlation  coefficient  can  be  either  positive  or  negative  depending  upon  whether  the 
atmosphere  is  unstable  or  stable,  respectively.  Generally,  it  is  accepted  to  have  an  absolute 
value  between  0.6  and  1.  Following  Claverie  et  al.  (Ref.  3),  L(W)WKD  uses  an  absolute 
value  of  0.8  as  the  exact  value  has  been  shown  to  have  little  influence  on  the  final  result. 
Claverie  et  al.  (Ref  3)  state  that  the  constant  of  proportionality,  Aj,,  can  vary  from  1  to  10"* 
in  the  case  of  strong  vapour  pressure  fluctuations.  Neglecting  the  strong  fluctuations,  they 
assume  it  has  a  value  of  1,  even  though  this  may  not  be  entirely  realistic.  Davidson  et  al. 
(Ref  24)  suggest  a  value  for  Aj,  of  0.6.  For  the  moment  L(W)WKD  will  follow  the 
recommendations  of  Claverie  et  al.  until  more  experimental  data  have  been  collected.  In  this 
respect,  A21  can  be  thought  of  as  another  experimental  parameter  to  be  determined. 
Consequently,  both  PIRAM  and  L(W)WKD  use  the  same  values  for  A21,  and  the  correlation 
coefficient. 


3.0  THE  LtWlWKD  PROGRAM 


Version  7.09  of  the  L(W)WKD  program  has  seen  a  number  of  important  changes  since 
Version  6. 1  was  discussed  in  DREV’s  earlier  report  (Ref  1).  The  principal  changes  include: 

1)  changes  to  the  modular  structure  and  their  functions  so  as  to  satisfy  the  Canada- 
France  agreement  on  a  common  set  of  fundamental  equations  and  options 

2)  the  capability  for  the  user  to  modify  certain  global  parameters  within  the 
program 

3)  the  capability  to  compile  the  program  in  three  different  modes  for  easier 
integration  into  modular  programs  such  as  IRBLEM  (Ref  8) 

and 

4)  the  ability  to  calculate  and  output  refractivity  profiles  for  multiple  sets  of 
meteorological  data  taken  at  differing  horizontal  ranges. 


Apart  from  these  four  changes,  which  are  discussed  in  more  detail  below,  many  other 
changes  have  been  made  since  Version  6.1.  These  changes  are  listed  in  the  C  file 
“LWWKD.C”,  of  which  its  relevant  section  is  given  for  reference  in  Appendix  B. 
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To  adapt  the  previous  L(W)WKD  program  so  as  to  integrate  the  common  set  of 
equations  and  options  agreed  upon  thru  an  exchange  between  Canada  and  France,  and  for 
other  reasons  as  well,  the  modular  structure  of  the  program,  and  many  of  its  functions  were 
modified.  Most  of  the  changes  relevant  to  the  Canada-France  agreement  were  made  to  the 
PROFILE.C  and  WEATHER.C  modules  and  their  header  files,  as  these  files  contain  all  of 
the  profile  and  meteorological  (weather)  functions. 

In  order  to  give  greater  flexibility  to  the  program,  the  program  now  looks  for  the 
initialization  file  “LWWKD.INF’.  If  it  doesn’t  exist  the  program  uses  default  values  for 
certain  global  parameters,  while  if  it  does  exist  the  program  will  use  the  user  defined 
parameters  contained  within  the  file.  This  addition  required  the  addition  of  the  module 
READINI.C  and  its  header  file.  Further  detail  about  the  use  of  the  initialization  file  are  give 
in  a  later  section. 

Earlier  versions  of  L(W)WKD  have  been  able  to  operate  in  BATCH  mode.  That  is,  the 
user  could  runlhe  program,  select  a  METDATA  file,  and  the  program  would  output  files 
containing  the  vertical  refractivity  profiles  for  a  desired  wavelength  band,  and  various 
vertical  meteorological  (temperature,  pressure,  etc.)  profiles.  However,  in  recent  years,  a 
variants  of  the  program  have  been  produced  so  that  they  could  be  integrated  as  modules  into 
modular  programs  such  as  IRBLEM.  As  a  result,  in  order  to  reduce  the  amount  of  effort 
required  to  update  these  versions  of  the  program,  the  same  seventeen  C  modules  and  eighteen 
header  files  that  currently  make  up  L(W)WKD,  can  be  compiled  into  three  different 
programs  by  simply  changing  the  INTERACT  and  WKDCN2  definitions  within  the 
LWWKD.H  header  file.  If  INTERACT  is  set  to  1  and  WKDCN2  set  to  0,  the  INTERAC¬ 
TIVE  version  of  the  program  LWWKD  is  created;  however,  when  INTERACT  is  set  to  0, 
two  different  BATCH  programs,  both  of  which  are  designed  to  be  integrated  into  IRBLEM 
(Ref.  8)  as  separate  modules,  are  created  depending  upon  whether  WKDCN2  is  set  to  0  or 
1.  If  WKDCN2  is  set  to  0,  a  BATCH  version  of  the  program  is  created  that  outputs  a  file 
containing  the  vertical  profile  of  the  refractive  index  and  its  gradients,  while  if  WKDCN2 
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is  set  to  1,  a  batch  version  of  the  program  is  created  that  outputs  a  file  containing  the  vertical 
profile  of  the  refractive  index  structure  parameter.  Both  BATCH  programs  operate  in  a 
similar  fashion  to  the  BATCH  mode  of  operation  within  the  INTERACTIVE  program,  in 
that  they  look  for  a  METDATA  file  (Ref  1),  processes  it,  and  create  appropriate  output  files 
for  use  in  modular  programs  such  as  IRBLEM.  More  details  about  the  three  compilation 
modes  are  given  in  a  later  section. 

The  ability  to  calculate  multiple  vertical  refractivity  profiles  using  meteorological  data 
at  differing  horizontal  ranges  is  not  an  addition  to  Version  7.09;  however,  it  wasn’t  present 
in  Version  6. 1 .  This  additional  feature  is  only  currently  available  with  the  INTERACTIVE 
version  of  the  program.  It  is  useful  when  the  vertical  refractivity  profiles  are  saved  as  an 
L(W)WKD  binary  data  file  and  that  file  is  used  with  a  compatible  ray-tracing  program  that 
can  accept  two-dimensional  (vertical  elevation  and  horizontal  range)  refractivity  data. 
Ref  25  shows  an  example  of  its  use  with  the  DREV  program  REFRACT. 

Version  7.09,  is  written  for  DOS  using  Version  3.0  of  Borland's  TURBO  C/C**.  The 
program  can  also  be  run  as  a  DOS  program  in  WINDOWS  or  OS/2.  Presently,  the  code  is 
composed  of  seventeen  C  modules  and  eighteen  user  defined  header  files  containing  ap¬ 
proximately  360  kb  of  code.  All  three  compilations  are  compiled  and  linked  using  the  large 
memory  model  and  produce  about  260  kb  of  executable  code.  The  names  of  the  modules  are 
LWWKD.C,  LWWKDBAT.C,  LWWKDINT.C  LWWKDSCR.C,  GRAPHMET.C, 
GRAPHREF.C,  GRAPHOTH.C,  PROFILE.C,  WEATHER.C,  READINI.C,  LBERROR.C, 
LFTLE.C,  L.GRAPH.C,  WIN.C,  MENU.C,  VALUEX.C,  and  MISC.C.  Similarly,  the 
header  files  afe  LWWKD.H,  LWWKDBAT.H,  LWWKDINT.H,  LWWKDSCR.H, 
GRAPHMET.H,  GRAPHREF.H,  GRAPHOTH.H,  PROFILE.H,  WEATHER.H, 
READINI.H,  LIBERROR.H,  LFILE.H,  L_GRAPH.H,  WIN.H,  MENU.H,  VALUEX.H, 
CDEF.H,  and  IBMKEYS.H.  LWWKD.C  is  the  entrance  module  and  contains  two  different 
main  routines  whose  use  depends  upon  the  setting  of  the  INTERACT  definition  within  the 
LWWKD.H  header  file.  Figures  3a,b  and  c  show  the  flow  of  the  L(W)WKD  program  in 
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both  INTERACTIVE  and  BATCH  compilation  modes. 

3.1  INTERACTIVE  Program 

If  the  INTERACTIVE  program  is  compiled,  the  user  can  interactively  control  and 
monitor  the  program’s  operation.  This  interactivity  is  achieved  through  the  use  of  menus 
whose  options  are  selected  using  the  keyboard,  and  windows  that  can  either  display 
information  to  the  user,  request  input  from  the  user,  or  allow  the  user  to  edit/modify 
information  used  by  the  program.  The  modules  WIN.C,  MENU.C  and  MISC.C  along  with 
the  header  files  IBMKEYS.H,  WIN.H,  MENU.H  and  CDEF.H  contain  all  the  routines  that 
manage  these  interactive  resources.  Further  information  regarding  these  resources  can  be 
found  in  Ref.  1 . 


3.1.1  LrWlWKD’s  Main  Screen 

Upon  executing  the  L(W)WKD  INTERACTIVE  program,  a  copyright  page  informs  the 
user  of  the  owner's  proprietary  rights,  the  program’ s  version  and  creation  date.  After  reading 
the  warning,  the  user  can  press  any  key  to  enter  the  program.  The  program  starts  by  looking 
for  the  LWWKD  initialization  file  “LWWKD.INr’,  and  reading  it  if  it  exists.  Table  n  shows 
the  general  format  for  the  initialization  file  and  the  default  values  for  all  the  user-definable 
parameters.  It  also  shows  the  minimum  and  maximum  permissible  values.  Any  improper 
value  causes  the  program  to  display  a  warning  message,  and  to  use  the  default  values  instead 
of  those  contained  within  the  initialization  file.  Subsequently,  it  initializes  other  default 
structures  within  the  program,  calculates  a  default  profile  parameter,  and  displays  the  Main 
screen  and  the  Program  menu. 

The  Main  screen  is  composed  of  two  text  windows  and  a  menu  window  (see  Fig.  4). 
The  top-left  "Weather  Measurements"  window  shows  the  default  weather  parameters  used 
by  L(W)WKD  to  calculate  the  profile  parameters  that  are  shown  in  the  top-right  "Weather 
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TABLE  n 

General  format  for  the  LWWKD.TNT  file 


Line# 

Variables 

Default 

Minimum 

Maximum 

CommentAExampIe 

1 

Comment 

-  all  comment  lines  must  be  at  the  beginning  of 
the  file,  and  they  must  start  with  a  #. 

-  for  example:  #  This  is  a  comment. 

2 

kel_cel 

273.15 

273 

274 

-  variable  to  go  between  Celsius  and  Kelvin 

-  for  example:  kel_cel  273. 1 

3 

Mair 

28.9 

28 

30 

-  variable  giving  the  mass  of  dry  air 

-  for  example:  Mair  28.8 

m 

Mwater 

18 

17 

20 

-  variable  giving  the  mass  of  water 

-  for  example:  Mwater  18.1 

5 

R 

8.3145 

8 

9 

-  variable  giving  the  Gas  Constant 

-  for  example:  R  8.31 

6 

g 

9.806 

9.5 

10 

-variable  giving  the  acceleration  due  to  gravity 
-  for  example:  g  9.8 

1 

gamm 

1.4 

1.35 

1.45 

-  variable  giving  the  ratio  of  specific  heat  at 
constant  pressure  to  that  at  constant  volume  for 
dry  air 

-for  example:  gamm  1 .4 1 

8 

cp 

1004 

1000 

1010 

-  variable  giving  the  specific  heat  at  constant 
pressure 

-  for  example:  cp  1004. 1 

9 

Kconst 

0.4 

0.35 

0.45 

-  variable  giving  von  Karman’s  constant 

-  for  example:  Kconst  0.41 

10 

p_adiab 

1000 

950 

1050 

-  variable  giving  the  adiabatic  reference 
pressure 

-  for  example:  p_adiab  1010.0 

11 

max_height 

50 

5 

100 

-  variable  giving  the  maximum  output  height 
for  vertical  profiles 

-  for  example:  max_height  20.0 

12 

pot_mode 

1 

0 

1 

-  variable  giving  the  potential  temperature 
mode  (0  -  use  dry  air  pressure;  1  -  use  total  air 
pressure) 

-  for  example:  pot_mode  0 

13 

scale_mode 

1 

0 

1 

-  variable  giving  the  scaling  function  mode  for 
stable  conditions  (0  -  LWKD;  1  -  KONDO) 

-  for  example:  scale_mode  0 

Note:  all  12  variables  do  not  have  to  be  listed  and  they  do  not  have  to  be  listed  in  the  above  order.  Thus 
the  user  only  has  to  include  those  variables  for  which  he  does  not  want  to  use  the  default  values. _ 
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Profiles"  window.  Three  changes  from  Ver.  6.1  can  be  noticed  on  the  Main  screen.  These 
are  the  inclusion  of  multiple  weather  parameter  data  sets  at  different  horizontal  ranges,  the 
addition  of  a  sea  surface  relative  humidity  to  each  data  set,  and  the  addition  of  the  ability  to 
choose  among  three  different  model  options.  Table  HI  shows  the  default  weather  parameters 
(not  tested  means  that  no  bounds  test  is  made).  Figure  4  also  shows  these  same  parameters 
along  with  the  resultant  vapour  pressure  for  the  chosen  data  set  as  shown  by  the  title 
“WEATHER  PARAMETERS  - 1"  on  the  left  side  of  the  screen.  At  the  bottom  of  the  figure, 
the  three  default  modes  for  the  options  Scaling  (Scaling  Equations),  Pot.  Temp.  (Potential 
Temperature),  and  Pressure  (Reference  Adiabatic  Pressure)  are  also  shown.  By  default,  the 
KONDO  version  of  the  scaling  equations  is  used  under  stable  conditions  (see  Eq.  3),  and  the 
total  atmospheric  pressure  (Tot.)  with  a  reference  pressure  of  1000  mb  are  used  in  Eq.  8  to 
determine  the  potential  temperature. 


TABLE m 

Parameter  Bounds  and  Defaults  in  Interactive  Mode 


Parameter/Units 

Minimum  Bound 

Maximum  Bound 

Default 

Horizontal  Range  (km) 

0 

100 

0 

Wave  Height  (m) 

0 

10 

0.25 

Wind  Sensor  Height  (m) 

0 

50 

15 

Wind  Speed  (m/s) 

0 

not  tested 

5 

T/P/RH  Sensor  Height  (m) 

0 

50 

15 

Air  Temperature  (°C) 

not  tested 

not  tested 

10 

Sea  Temperature  (°C) 

not  tested 

not  tested 

5 

Atmospheric  Pressure  (mb) 

0 

not  tested 

1010 

Relative  Humidity 

0 

1 

0.8 

Sea  Surface  Relative  Humidity  (%) 

0 

1 

0.982 

The  L(W)WKD  profile  parameters  calculated  using  these  initial  values  are  displayed 
in  the  right-hand  side  window  of  Fig.  4.  It  shows  the  calculated  roughness  heights  and 
scaling  constants  used  for  the  wind  speed,  virtual  potential  temperature  and  the  specific 
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humidity  profiles.  It  also  shows  the  wind  speed,  air  temperature,  relative  humidity,  and  water 
vapour  pressure  calculated  using  these  profiles,  along  with  the  Monin-Obukhov  length  and 
the  number  of  iterations  that  were  required  for  L(W)WKD  to  obtain  the  solution. 

The  Program  menu  (Fig.  3a)  is  the  point  of  entry  into  the  L(W)WKD  program  and  is 
composed  of  three  function  keys  F 1  -Interactive,  F2-Batch  and  F5-EXIT.  Pressing  the  F5  key 
opens  a  dialog  window  that  requires  a  user  response  to  the  following  question:  "Exit 
program  (Y/N)?".  A  No  response  leaves  the  user  inside  the  program  while  a  Yes  response 
clears  the  screen,  closes  the  program,  and  returns  the  user  to  the  operating  system.  Using  the 
FI -Interactive  function  key  opens  a  message  window  that  displays  “Use  the  Arrow  keys  to 
cycle  between  data  sets”  for  several  seconds  before  the  user  is  presented  with  the  Main  menu 
and  finds  himself  in  the  Interactive  mode  of  operation.  Likewise,  pressing  the  F2-Batch 
function  puts  the  user  in  the  Batch  mode  of  operation.  The  details  of  using  either  operational 
mode  are  fully  described  in  the  following  two  sections. 

3.1.2  Interactive  Mode 

After  the  FI -Interactive  function  key  of  the  Program  menu  has  been  selected,  another 
window  is  added  to  the  screen  and  the  Program  menu  is  replaced  by  the  Main  menu  (see 
Fig.  5).  The  new,  initially  empty  window  is  the  "Other  Weather  Measurements"  window. 
The  purpose  of  the  window  is  to  display  any  other  weather  measurements  that  the  user 
wishes  to  have  included  in  any  of  the  meteorological  graphs.  More  details  related  to  this 
window,  along  with  the  actions  initiated  by  the  four  function  keys;  Fl-Modify,  F2-Load 
Data,  F3-Save  Data,  and  F4-Output  Data,  of  the  Main  menu  are  given  in  the  following 
subsections.  The  F5-RETURN  key  returns  the  user  to  the  Program  menu. 
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3.1.2.1  FI -Modify 

The  FI -Modify  function  key  allows  the  user  to  modify  the  measured  parameters  so  as 
to  produce  new  profile  parameters,  to  modify  the  three  options,  and  to  add/delete/modify 
values  displayed  in  the  "Other  Weather  Measurements"  window.  Upon  pushing  the  key  a 
new  menu,  the  Modify  menu,  appears.  It  is  composed  of  the  four  function  keys  Fl-Weather, 
F2-Options,  F3-Other;Wind,  F4-Other;Temp  and  F5-RETURN  (see  Fig.  3b).  Pressing  the 
F5-RETURN  key  returns  the  user  to  the  Main  menu.  The  Fl-Weather  key  allows  the  user 
to  modify  the  weather  parameters,  the  F2  key  allows  the  user  to  modify  the  various  options, 
while  the  F3  and  F4  keys  allow  the  user  to  modify  the  wind  parameters  or  the  tempera¬ 
ture/relative  humidity  parameters  belonging  to  the  "Other  Weather  Measurements"  window. 

3. 1.2.2  Fl-Weather 

Pressing  the  Fl-Weather  key  presents  the  user  with  the  Weather  menu  (see  Fig.  3b). 
It  is  composed  of  the  following  four  function  keys:  FI -Edit,  F2-Add,  F3-Delete,  F4-RESET, 
F5-ABORT,  and  F6-RETURN.  The  F2-Add  key  is  only  displayed  while  the  number  of  data 
sets  is  less  than  10,  and  the  F3-Delete  key  is  only  displayed  when  there  is  more  than  1  data 
set.  The  F4-RESET  and  F5-ABORT  keys  are  only  present  after  a  weather  set  has  been 
edited.  Pressing  the  F6-RETURN  key  opens  the  following  dialog  window:  "Return  to  Main 
Menu?  (Y/N):".  A  No  response  leaves  the  user  at  the  same  place  within  the  program,  while 
a  Yes  response  returns  the  user  to  the  Modify  menu. 

Pushing  the  FI -Edit  key  opens  an  edit  window  through  which  the  user  can  modify  the 
horizontal  range  for  the  measurements,  the  Hj/3  wave  height,  the  wind  speed  and  the  wind 
sensor  height,  the  sea  temperature,  sea  surface  relative  humidity,  and  the  air  pressure,  air 
temperature,  relative  humidity,  and  their  sensor  height.  After  the  user  has  accepted  the 
changes  to  these  weather  parameters,  the  program  verifies  that  they  are  within  the  accepted 
bounds  given  by  Table  HI.  Any  parameters  outside  these  limits  are  reset  to  their  previously 
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held  values,  new  profile  parameters  are  calculated,  and  the  new  weather  and  profile 
parameters  are  displayed  on  the  main  screen.  On  the  other  hand,  if  the  profile  parameters 
cannot  be  calculated,  the  following  message,  "Could  not  fit  new  weather  parameters!!",  is 
displayed  in  a  message  window  for  two  seconds,  and  the  newly  edited  weather  parameters 
are  reset  to  their  previous  values. 

After  accepting  the  modifications  and  returning  to  the  Weather  menu,  the  F4-RESET 
and  F5- ABORT  keys  are  available.  The  F4-RESET  key  allows  the  user  to  reset  the  weather 
parameters  to  the  values  they  had  before  they  were  last  edited,  while  the  F5-ABORT  key 
opens  the  dialog  window:  “Reset  parameters  and  QUIT?  (Y/N):”.  A  No  response  leaves  the 
user  at  the  same  place  within  the  program,  while  a  Yes  response  resets  the  weather 
parameters  and  returns  the  user  to  the  Modify  menu  (i.e.,  it  does  a  RESET  and  a  RETURN). 

Pushing  the  F2-Add  key,  when  it  is  available,  displays  the  message,  "New  weather 
parameter  set  added!!",  for  two  seconds.  The  new  data  set  is  created  using  the  default 
weather  parameters  at  a  horizontal  range  of  0  km.  To  edit  or  delete  the  new  data  set,  use  the 
left/right  arrow  keys  to  select  the  correct  data  set  and  push  F 1  -Edit  or  F3-Delete,  respectively. 
The  operator  has  to  select  the  proper  weather  parameter  set  as  the  program  automatically 
reorders  the  sets  by  increasing  horizontal  range. 

Pushing  the  F3-Delete  key,  when  it  is  available,  displays  the  dialog  window:  “Delete 
current  weather  parameter  set?  (Y/N):”.  If  the  correct  set  has  not  been  chosen,  make  a  No 
response,  choose  the  correct  data  set  using  the  left/right  arrow  keys,  and  push  the  F3-Delete 
key  again.  Choosing  the  Yes  response  automatically  deletes  the  current  data  set,  reorders  the 
data  sets,  and  displays  one  of  the  remaining  data  sets  on  the  screen.  This  process  can  not  be 
undone. 
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3.1.2.3  F2-Options 

Pressing  this  new  function  key  presents  the  user  with  an  edit  window  that  allows  the 
user  to  modify  the  three  modes  for  the  Scaling  Equations,  the  Potential  Temperature,  and 
the  Reference  Pressure.  By  default  the  Scaling  Equations  parameter  is  set  to  1,  which 
corresponds  to  the  Kondo  version  of  the  gradient  (scaling)  equations  used  under  stable 
conditions  (see  Eq.  3),  the  Potential  Temperature  parameter  is  set  to  1,  which  corresponds 
to  using  the  total  atmospheric  pressure  in  Eq.  8,  and  the  Reference  Pressure  is  set  to 
1(X)0  mb  for  use  in  Eq.  8.  Changing  the  Scaling  Equations  parameter  to  0  allows  the  user 
to  use  the  WKD  version  under  stable  conditions,  changing  the  Potential  Temperature 
parameter  to  0  allows  the  user  to  use  the  dry-air  atmospheric  pressure  in  Eq.  8  instead  of  the 
total  air  pressure,  and  the  Reference  Pressure  can  take  any  value  greater  than  0  for  use  in 
Eq.  8.  Any  parameters  outside  these  limits  are  reset  to  their  previously  held  values,  new 
profile  parameters  are  calculated,  and  the  new  weather  and  profile  parameters  are  displayed 
on  the  main  screen.  On  the  other  hand,  if  the  profile  parameters  cannot  be  calculated,  the 
following  message,  "Could  not  fit  using  new  options ! ! ",  is  displayed  in  a  message  window 
for  two  seconds  and  the  newly  edited  option  parameters  are  reset  to  their  previous  values. 

After  accepting  the  modifications  and  returning  to  the  Options  menu,  the  F2-RESET 
and  F3-ABORT  keys  are  available.  The  F2-RESET  key  allows  the  user  to  reset  the  option 
parameters  to  the  values  they  had  before  they  were  last  edited  while  the  F3-ABORT  key 
opens  the  dialog  window:  "Reset  options  and  QUIT?  (Y/N); A  No  response  leaves  the  user 
at  the  same  place  within  the  program  while  a  Yes  response  resets  the  weather  parameters  and 
returns  the  user  to  the  Modify  menu  (i.e.,  it  does  a  RESET  and  a  RETURN). 

3. 1.2.4  F3-Other:Wind  &  F4-Other:Meteo 


Selecting  either  of  these  two  function  keys  presents  the  user  with  either  the  Wind  menu 
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or  the  Meteo  menu  (see  Fig.  3b).  Both  menus  are  identical  and  are  composed  of  the 
following  four  function  keys:  FI -Edit,  F2-Add,  F3-Delete  and  F4-RETURN.  However,  the 
FI,  F2  and  F3  keys  are  not  always  displayed.  Pressing  the  F4-RETURN  key  opens  the 
following  dialog  window:  "Return  to  Main  Menu?  (Y/N):".  A  No  response  leaves  the  user 
at  the  same  place  within  the  program  while  a  Yes  response  returns  the  user  to  the  Modify 
menu.  The  functions  of  the  other  keys  are  given  in  the  following  subsections. 

3.1.2.4.1  Fl-Edit 

The  F 1  -Edit  key  is  only  displayed  when  there  is  data  that  can  be  edited.  Pushing  the  key 
opens  an  edit  window  which  requests  the  user  to  enter  the  number  (1  to  8)  associated  with 
the  data  set  to  be  modified.  An  illegal  entry  produces  the  following  two  second  message 
window,  "Input  out  of  range! ! ",  and  then  allows  the  user  to  enter  a  new  number.  Once  the 
data  set  has  been  selected,  an  edit  window  is  opened  and  the  user  can  modify  the  associated 
data.  If  the  Other;Wind  function  key  had  previously  been  selected,  the  user  can  modify  the 
wind  speed  and  the  wind  sensor  height,  while  if  the  Other;Meteo  function  key  had  previously 
been  selected,  the  user  can  modify  the  air  temperature,  relative  humidity,  and  the  sensor 
height.  After  the  user  has  accepted  the  changes  to  the  edited  parameters,  the  program  verifies 
that  they  are  within  the  same  bounds  given  by  Table  HI,  before  accepting  them  and 
displaying  them  on  the  main  screen.  The  data  sets  are  always  displayed  in  increasing  order 
with  respect  to  their  sensor  heights. 


3.1. 2.4.2  F2-Add 

The  F2-Add  key  is  only  displayed  when  there  is  still  space  for  another  data  set  (the 
maximum  number  is  8).  If  there  is,  pushing  the  key  opens  an  edit  window  through  which 
the  user  can  modify  the  parameters  belonging  to  the  appropriate  default  data  set.  If  the 
Other;Wind  function  key  had  previously  been  selected,  the  user  can  add  a  new  wind  speed 
and  wind  sensor  height,  while  if  the  Other;Meteo  function  key  had  previously  been  selected. 
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the  user  can  add  a  new  air  temperature,  relative  humidity  and  meteo  sensor  height.  After  the 
user  has  accepted  the  new  data  set,  the  program  verifies  that  they  are  within  the  limits  given 
by  Table  IH,  and  displays  them  on  the  main  screen.  Again,  the  data  sets  are  listed  in 
increasing  ordered  with  respect  to  the  sensor  heights. 

3.1.2.4.3  F3-Delete 

The  F3-Delete  key  is  only  displayed  when  there  is  data  that  can  be  deleted.  Pushing  the 
key  opens  an  edit  window  which  requests  the  user  to  enter  a  number,  from  1  to  8,  of  the  data 
set  to  be  deleted.  An  illegal  entry  produces  the  following  two  second  message  window, 
"Input  out  of  range! !",  and  then  allows  the  user  to  enter  a  new  number.  Once  the  data  set  is 
selected,  the  following  dialog  window  is  displayed:  "Delete  data  set  (Y/N)?".  A  No  response 
aborts  the  delete  and  returns  the  user  to  the  preceding  menu  (Wind  or  Meteo)  while  a  Yes 
response  deletes  the  data  set  before  returning  the  user  to  the  preceding  menu.  The  data  sets 
are  then  reordered  with  respect  to  the  sensor  heights  before  they  are  redisplayed  on  the  main 
screen. 

3. 1.2.5  F3-Save  Data 

The  F3-Save  Data  function  key  allows  the  user  to  save  certain  information  to  a  special 
L(W)WKD  binary  file  having  the  extension  "Iwd".  The  file’s  most  recent  structure  is  shown 
in  Table  IV  where  the  first  column  gives  the  data  type  (integer,  float,  etc.),  the  second  gives 
the  number  of  bytes  it  requires,  and  the  third  column  describes  the  data  type’s  significance. 
This  file  is  a  general  purpose  file  that  includes  each  set  of  measured  weather  parameters,  and 
their  associated  profile  parameters  and  refractivity  profiles.  The  refract! vity  profiles  are 
designed  for  use  as  input  into  ray-tracing  algorithms.  The  file  can  be  extended  to  contain 
other  data  while  remaining  backward  compatible  by  changing  the  file  identification  number. 
In  this  way  each  file  identification  number  can  be  uniquely  associated  with  a  specific 
L(W)WKD  version’ s  file  format  such  that  any  program  using  such  files  always  know  exactly 


UNCLASSIFIED 

41 


how  to  retrieve  the  required  data. 

Upon  pressing  the  key,  the  user  is  presented  with  the  Wave  Selection  window  (see 
Fig.  3a).  It  is  a  selection  window  from  which  the  user  can  choose  from  one  of  five 
wavelength  bands  (Visible,  Near-IR,  Mid-IR,  Far-IR,  or  Radar).  After  choosing  the 
wavelength  band,  the  user  is  presented  with  the  Model  Selection  window.  This  is  a 
selection  window  from  which  the  user  must  choose  to  incorporate,  or  not,  the  Wavy  Model 
for  calculating  the  refractivity  profiles.  Once  both  selections  are  made,  a  Save  File  window 
is  produced  showing  a  listing  of  all  the  files  in  the  current  directory  having  the  extension 
"Iwd".  Scroll  through  this  listing  until  an  edit  window  appears  and  enter  a  filename.  Do  not 
include  the  extension,  as  the  extension  "Iwd"  is  automatically  appended.  Once  the  filename 
has  been  accepted,  the  user  is  presented  with  a  second  edit  window.  This  window  allows  the 
user  to  create  a  unique  file  header  for  every  L(W)WKD  file.  After  creating  the  file  header, 
the  file  is  created,  saved  to  the  current  directory,  and  the  program  displays  the  message: 
"Writing  to  file . Each  file  requires  at  least  4.5  kb. 


TABLE  IV 

Format  for  the  LfWlWKD  binary  file 


Bytes 

Comments 

2 

File  identification  number  =  4.  Previous  versions  used  0, 1,2,  &  3. 

2 

#  of  bytes  in  the  file  header  including  the  null  character  =  N. 

3.  String 

N 

N  characters  belonging  to  the  file  header. 

2 

#  of  bytes  in  the  wavelength  band’s  name  including  the  null  character  =  M. 

5.  String 

M 

M  characters  belonging  to  the  wavelength  band’s  name 

-  VISIBLE,  NEAR-INFRARED,  MID-B'IFRARED,  FAR-INFRARED  or 

RADAR. 

2 

#  of  weather  sets  (maximum  is  10) 

7.  Character 

1 

Scaling  Equation  Option  (default  is  1) 

8.  Character 

1 

Potential  Temperature  Option  (default  is  1) 

9.  Double 

8 

Reference  Pressure  Option  (default  is  1000  mb) 

The  following  4  outputs  (10  to  13)  are  written  for  each  weather  set. 
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2 

#  of  bytes  in  the  weather  parameter  structure  (PARAM). 

ll.PARAM 
(11  floats) 

32 

PARAM  contains  the  following  1 1  weather  parameters 
-  horizontal  range,  H,;3  wave  height,  wind  sensor  height,  wind  speed,  meteo 
sensor  height,  air  temperature,  sea  tenperature,  air  pressure,  relative  humid¬ 
ity,  sea  surface  rel.  humidity,  and  the  sea  surface  wind  speed  (=  0). 

2 

#  of  bytes  in  the  weather  profile  structure  (W_STRUCT). 

13.  W_STRUCT 
(1  unsigned) 
(10  double) 

81 

W_STRUCT  contains  the  following  profile  parameters 
-  number  of  iterations  (unsigned  int),  Monin-Obukhov  length,  wind  speed 
scale  and  roughness  height,  wind  speed,  virtual  potential  temperature  scale 
and  roughness  height,  air  temperature,  water  vapour  density  scale  and 
roughness  height,  and  relative  humidity. 

14.  Float 

4 

Minimum  height  =  H,„/2. 

15.  Float 

4 

Maximum  height  =  50. 

2 

#  of  bytes  in  the  Other  Wind  weather  structure  (W_PARAM). 

2 

#  of  Other  Wind  weather  structure  data  sets. 

Output  18  is  done  for  each  Other  Wind  data  set. 

8 

W  PARAM  contains  a  wind  sensor  height  and  a  wind  speed. 

2 

#  of  bytes  in  the  Other  Meteo  weather  structure  (M_PARAM). 

2 

#  of  Other  Meteo  weather  structure  data  sets. 

Output  21  is  done  for  each  Other  Wind  data  set. 

21.M_PARAM 

12 

M_PARAM  contains  a  meteo  sensor  height,  an  air  temperature  and  a 
relative  humidity. 

The  following  output  is  written  for  each  weather  set. 

22.  Data  Record 
(2  float) 

(3  double) 

i 

32 

Each  data  record  contains  the  range  (L  =  0  km),  the  height  (m)  above  the 
MWL  (h),  the  refractivity  (Z),  and  two  derivatives  (dZ/dh  and  dZ/dL(=0)) 
at  the  height  h.  The  height  (h)  is  calculated  iteratively  as  follows:  h  =  h  + 
0.05  +  0.025*h  and  records  are  written  until  h  =  max.  height.  The  refrac- 
tivities  are  calculated  for  the  chosen  wave  band. 

3. 1.2.6  F2-Load  Data 


The  F2-Load  Data  function  key  allows  the  user  to  load  the  meteorological  parameters 
previously  saved  in  an  L(W)WKD  binary  file  from  the  current  directory.  When  the  key  is 
pressed,  the  user  is  presented  with  a  Load  File  window  showing  a  numbered  list  of  all  the 
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files  in  the  current  directory  with  the  extension  "Iwd".  Make  a  note  of  the  number  associated 
with  the  file  to  be  loaded  while  scrolling  through  the  listing.  When  the  edit  window  appears, 
enter  the  number.  If  the  number  is  valid,  the  file  is  loaded,  and  displays  the  file’ s  header  and 
the  wavelength  band  for  two  seconds  each  before  displaying  the  message  "File  Loaded....". 
It  should  be  noted,  that  all  previous  binary  file  formats  can  be  read,  and  that  the  refractivity 
profile  data  records  are  neither  read  nor  required  by  the  program. 


3. 1.2.7  F4-OutDut  Data 


Pressing  the  F4-Output  Data  key  of  the  Main  menu  presents  the  user  with  the  Output 
menu  (see  Fig.  3c).  This  menu  consists  of  the  five  function  keys:  Fl-Meteo  Screen,  F2- 
Meteo  Profile,  F3-Ref  Profile,  F4-Other  Graphs,  and  F5-RETURN.  Pressing  the  F5- 
RETURN  key  returns  the  user  to  the  Main  menu  while  use  of  the  other  keys  allows  the  user 
to  produce  various  textual  or  graphical  outputs.  These  outputs  are  discussed  in  the  following 
subsections. 


3. 1.2.7. 1  Fl-Meteo  Screen 

Choosing  the  Fl-Meteo  Screen  key  allows  the  user  to  output  the  data  displayed  on  the 
Main  screen  to  a  file  in  either  ASCII  or  PostScript  format.  These  files  can  then  be  sent  to 
a  printer  to  produce  a  hardcopy  of  the  Main  screen  or  as  input  for  other  software.  Upon 
pressing  the  key,  the  user  is  presented  with  a  File  Selection  window.  This  selection  window 
contains  two  choices:  PostScript  or  ASCII  File.  After  selecting  the  file  type,  a  Save  File 
window  is  produced  showing  a  listing  of  all  the  files  in  the  current  directory  with  the 
appropriate  extension  ("psc"  or  "asc").  Scroll  through  the  list  until  the  edit  window  appears 
and  enter  a  filename  without  an  extension.  Once  the  filename  is  accepted,  the  file  is  created, 

saved  to  the  current  directory,  and  the  program  displays  the  message:  "File  Written . ", 

before  returning  to  the  File  Selection  window.  The  PostScript  file  requires  5.5  kb,  and  the 
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ASCn  file  1.5  kb. 


3.1. 2.7.2  F2-Meteo  Profile 

Choosing  the  F2-Meteo  Profile  key  allows  the  user  to  graphically  output  various  vertical 
profiles,  for  six  different  meteorological  parameters,  to  the  monitor  or  a  PostScript  file. 
Upon  pressing  the  key,  the  user  is  presented  with  the  Meteo  Selection  window.  This 
selection  window  has  the  following  options:  Wind  Speed,  Air  Temp.,  Virt.  Pot.  Temp.,  Total 
Pressure,  Rel.  Humidity,  Vapour  Pressure,  Spec.  Humidity,  Vapour  Density,  and  ALL- 
ASCII.  The  ALL-ASCn  option  allows  the  user  to  create  an  ASCII  file  containing  the 
profiles  for  all  eight  parameters.  The  file’s  general  format  is  given  in  Table  V.  After 
selecting  any  of  the  above  choices,  the  user  is  presented  with  the  Model  Selection  window. 
Once  the  model  (Non-Wavy  or  Wavy)  for  the  profile  calculation  is  chosen,  the  Height  Edit 
window  appears.  It  is  an  edit  window  requiring  the  user  to  enter  the  maximum  height  to 
which  the  profiles  are  to  be  calculated.  The  default  value  is  given  by  max_height  as  given 
by  Table  H.  After  entering  this  parameter  the  program  follows  one  of  two  paths. 


TABLE  V 

General  format  for  all  ASCII  files 


Line  # 

Variable  (Format) 

Comment 

1. 

File  identification  number  (integer) 

Identifies  the  file’s  origins  (=  10) 

2. 

File  description  (string) 

States  the  type  of  file  and  L(W)WKD  version 

3. 

Waveband  descriptor  (string) 

This  line  is  optional.  Not  all  files  use  it. 

4. 

Min.  and  max.  height  (reals) 

This  line  is  optional.  Not  all  files  use  it. 

5. 

Fixed  profile  height  ( real) 

This  line  is  optional.  Not  all  files  use  it. 

6. 

Headers  for  each  column  (string) 

The  headers  describe  the  structure  of  each  data  record 

■■ 

Underline  for  each  header  (string) 

8. 

Data  record  (all  real) 

Currently,  the  largest  record  contains  9  consecutive  real 
parameters 
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If  the  ALL-ASCn  option  was  chosen,  a  Save  File  window  is  displayed  showing  a  listing 
of  all  the  files  in  the  current  directory  with  the  extension  "asc".  Scroll  through  the  listing 
until  the  edit  window  appears,  and  enter  a  filename  without  an  extension.  After  the  filename 
is  accepted,  the  file  is  created,  saved  to  the  current  directory,  and  the  program  displays  the 

message:  "File  Written . ",  before  returning  to  the  Meteo  Selection  window.  The  ASCII 

file  requires  about  9  kb  for  the  default  height  of  50  m  with  a  step  size  of  0.5  m. 

If  any  of  the  other  options  were  chosen,  the  user  is  shown  the  Display  Selection 
window.  It  is  a  selection  window  with  the  following  two  display  options:  Monitor  or 
PostScript.  Choosing  Monitor  will  create  a  graphical  display  on  the  screen.  To  remove  the 
graph  from  the  screen  and  be  returned  to  the  Meteo  Selection  window,  the  user  may  strike 
any  key  on  the  keyboard.  If  PostScript  is  chosen,  a  Save  file  window  is  produced  showing 
a  listing  of  all  the  files  in  the  current  directory  with  the  extension  "psc".  Scroll  through  the 
list  until  the  edit  window  appears  and  enter  a  filename  without  an  extension.  Once  the 
filename  is  accepted,  the  graph  is  written  to  the  file  and  saved  to  the  current  directory,  before 
returning  the  user  to  the  Meteo  Selection  window.  The  PostScript  file  requires  close  to  5  kb 
under  the  default  conditions.  Any  weather  data  present  in  the  "Other  Weather  Measure¬ 
ments"  window  is  also  plotted  on  these  graphs;  however,  they  are  not  included  in  the  ASCII 
files. 


3.1.2.7.3  F3-Index  Profile 

Choosing  the  F2-Ref.  Profile  key  allows  the  user  to  graphically  output  XY  plots  of  either 
the  refractivity  (Z),  the  derivative  of  the  refractivity  with  height  (dZ/dh),  or  the  second 
derivative  of  the  refractivity  with  height  (d^Z/dh^)  for  any  of  the  wavelength  bands  to  the 
monitor,  or  to  a  PostScript  file  for  different  X-parameters.  The  possible  X-parameters  are 
the  height  above  the  MWL,  the  ASTD  (sea  temperature  fixed),  the  wind  speed,  and  the 
relative  humidity.  All  plots  can  also  be  done  with  or  without  the  Wavy  model. 
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Upon  pressing  the  key,  the  user  is  presented  with  the  Refractivity  Y  Selection  window. 
This  is  a  selection  window  having  the  following  options:  Refractivity,  d(Ref.)/d(H), 
dd(Ref.)/dd(H),  and  ALL- ASCII.  The  ALL- ASCII  option  allows  the  user  to  create  an  ASCII 
file  containing  the  values  for  all  three  refractivity  profiles  using  the  same  general  format 
shown  in  Table  V.  After  pressing  the  key,  the  user  is  presented  with  the  Refractivity  X 
Selection  window.  It’s  a  selection  window  with  the  following  options:  Height,  ASTD  (sea 
fixed),  Wind  Speed,  and  Rel.  Humidity.  Once  it  has  been  selected,  the  Wave  Selection 
window,  and  then  the  Model  Selection  window  are  displayed  to  the  user.  After  choosing 
both  the  wave  band  and  whether  or  not  to  use  the  Wavy  model,  the  Height  Edit  window 
appears.  If  the  X-dependent  parameter  is  "Height",  enter  the  maximum  height  (< 
max_height  given  by  Table  II)  to  which  the  appropriate  Y-parameter  profile  is  to  be  plotted. 
Otherwise,  enter  the  height  (<  max_height)  at  which  the  appropriate  Y-parameter  is  to  be 
calculated  for  the  chosen  X-dependent  parameter.  After  choosing  this  height,  the  program 
follows  one  of  the  two  following  paths. 

If  the  ALL-ASCn  option  was  chosen,  a  Save  File  window  is  displayed  showing  a  listing 
of  all  the  files  in  the  current  directory  with  the  extension  "asc".  Scroll  through  the  list  until 
the  edit  window  appears,  and  enter  a  filename  without  an  extension.  After  the  filename  is 
accepted,  the  file  is  created,  saved  to  the  current  directory,  and  the  program  displays  the 

message:  "File  Written . ",  before  returning  to  the  Ref.  Y  Selection  window.  The  ASCII 

file  usually  requires  less  them  10  kb. 

If  any  other  option  was  chosen,  the  user  is  shown  the  Display  Selection  window. 
Choosing  Monitor  creates  a  graphical  display  on  the  screen.  To  remove  it  from  the  screen 
and  be  returned  to  the  Wave  Selection  window,  the  user  may  strike  any  key  on  the  keyboard. 
If  PostScript  is  chosen,  a  Save  File  window  is  produced  showing  a  listing  of  all  the  files  in 
the  current  directory  with  the  extension  "psc".  Scroll  through  the  list  until  the  edit  window 
appears,  and  then  enter  a  filename  without  an  extension.  Once  the  filename  is  accepted,  the 
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graph  is  written  to  the  file  and  saved  to  the  current  directory  before  returning  the  user  to  the 
Ref.  Y  Selection  window.  The  PostScript  file  usually  requires  less  than  6  kb. 

3.1. 2.7.4  F4-Qther  Graphs 

Choosing  the  F4-Other  Graphs  key  allows  the  user  to  graphically  output  various  XY 
plots  to  the  monitor,  or  to  a  PostScript  file  using  either  the  Non-Wavy  or  Wavy  model. 
ASCn  files,  using  the  same  general  format  given  in  Table  V,  can  also  be  created.  Upon 
pressing  the  key,  the  user  is  presented  with  the  Graphic  Selection  window  (see  Fig.  3c).  It’s 
a  selection  window  with  the  following  two  general  options  for  the  Y-axis:  Structure 
Parameters  and  Profile  Parameters. 

Choosing  Structure  Parameters  allows  the  user  to  create  XY  plots  for  four  structure 
parameters  as  functions  of  six  different  parameters  for  any  of  the  five  wavelength  bands.  The 
possible  structure  parameters  are  the  refractivity  stmcture  parameter,  the  temperature 
structure  parameter,  the  vapour  density  structure  parameter,  and  the  temperature-vapour 
density  structure  parameter.  The  possible  dependent  X-axis  parameters  are  the  height  above 
the  MWL,  the  ASTD  (sea  temperature  fixed),  the  ASTD  (air  temperature  fixed),  the  wind 
speed,  the  relative  humidity,  and  the  H1/3  wave  height. 

Choosing  Profile  Parameters  allows  the  user  to  create  XY  plots  of  the  various  weather 
profile  parameters  versus  the  ASTD  (sea  temperature  fixed),  the  ASTD  (air  temperature 
fixed),  the  wind  speed,  the  relative  humidity,  the  reference  pressure  (adiabatic  pressure),  or 
von  Karmen’s  constant.  The  possible  weather  profiles  are  the  Monin-Obukhov  length,  the 
wind  scale  parameter,  the  wind  roughness  height,  the  temperature  scale  parameter,  the 
temperature  roughness  height,  the  humidity  scale  parameter,  and  the  humidity  roughness 
height.  As  the  process  involved  in  creating  these  graphs  or  files  are  nearly  identical  to  that 
for  creating  graphs  of  the  refractivity  profiles,  only  their  differences  will  be  described. 
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If  the  user  chooses  to  create  graphs  of  a  structure  parameter,  the  user  will  discover  that 
the  process  is  nearly  identical  to  that  for  creating  graphs  of  the  refractivity  profiles.  The 
differences  are  only  with  respect  to  the  available  choices  for  the  Y-axis  and  X-axis 
parameters.  On  the  other  hand,  if  the  user  chooses  to  create  graphs  of  a  profile  parameter, 
the  user  will  notice  the  following  differences.  After  having  entered  a  choice  for  the  Y- 
parameter  and  the  X-parameter,  the  user  will  not  be  asked  to  select  a  wave  band  or  whether 
or  not  to  use  the  Wavy  model.  This  is  because  these  parameters  are  independent  of  these  two 
options.  After  displaying  the  selected  XY  plot  on  the  monitor,  saving  the  selected  XY  plot 
to  a  PostScript  file,  or  saving  the  entire  set  to  an  ASCII  file,  the  program  returns  to  the 
appropriate  Y  Selection  window.  The  PostScript  files  use  less  than  5  kb  while  the  ASCII 
files  use  less  than  10  kb. 


3.1.3  Batch  Mode 

Pressing  the  F2-Batch  function  key  of  the  Program  menu  allows  the  user  to  produce 
weather,  structure  parameter,  and  refractivity  output  files  from  an  input  file  containing  up  to 
999  sets  (records)  of  weather  measurements.  Upon  selecting  this  mode,  three  new  windows 
are  added  to  the  bottom  of  the  Main  screen.  These  are  the  "Warnings",  "Status"  and 
“METDATA  File”  windows  (see  Fig.  6).  The  "Warnings"  window  informs  the  user  of  any 
problems  associated  with  analyzing  a  specific  weather  record  while  the  "Status"  window 
shows  the  user  which  weather  record  is  currently  being  processed.  The  “METDATA  File” 
window  is  an  edit  window  containing  the  name  of  the  default  input  ASCII  file: 
"METDATA".  The  user  may  change  the  name  or  accept  it.  The  format  of  this  file  is  given 
in  Table  VI.  If  a  valid  METDATA  file  is  not  entered,  the  user  is  returned  to  the  Program 
menu;  otherwise,  the  user  is  presented  with  the  Wave  Selection  window,  and  then  the  Model 
Selection  window.  After  having  selected  both  the  wavelength  band,  and  whether  or  not  to 
incorporate  the  Wavy  model,  the  Processing  menu  (see  Fig.  3a)  is  displayed  and  the  program 
begins  to  process  the  weather  records  belonging  to  the  METDATA  file.  The  Processing 
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TABLE  VT 

General  format  for  a  METDATA  ASCII  file  -  a  space  separates^eadijf^^ 


Line# 

Variables 

Comment/Example 

1 

String  comment  line 

-each  comment  line  must  start  with  a  # 

-  for  example:  #This  is  a  comment. 

2 

5  integers 

-day,  month,  hour,  minute,  and  measurement  interval  (min.) 

-  for  example:  20  2  8  20  5  (they  are  often  set  to  zero) 

3 

Integer  (N„) 

-number  of  weather  conditions  in  the  file 
-  for  example:  10 

4... 

Read  weather  condition 
(read  1  line  at  a  time) 

-each  line  contains  22  variables  in  the  following  order:  air 
temperature  measurement  height  (m),  atmospheric  pressure 
(mbar),  air  temperature  (“  C),  relative  humidity  (%),  solar 
radiation*,  wind  speed,  wind  direction*,  water  temperature 
(°C),  wind  speed  measurement  height  (m).  Hi, 3  wave  height 
(m),  visibility  (km)*,  humidity  measurement  height  (m)*, 
neutral  drag  coefficient  (A)*,  neutral  drag  coefficient  (B)*, 
neutral  Stanton  number*,  neutral  Dalton  number*,  wavelength 
(pm)*,  pressure  measurement  height  (m)*,  air  mass  parame¬ 
ter*,  and  24hr  wind  speed*,  precipitation  rate  (mm/hr)*,  and 
precipitation  type  (0-4)*. 

*  parameter  is  not  required  by  L(W)WKD 
-  for  example 

15.0  1010.0  10.0  80.0  500.0  5.0  270  12.0  20.0  1.5  9.0, 0, 0, 0, 

0, 0, 0, 0, 10,3.0  0  0 

menu  contains  two  active  function  keys:  El -Pause/Restart  and  F5-STOP.  Pressing  the  FI 
function  key,  which  initially  displays  "Pause",  causes  the  data  processing  to  be  stopped, 
changes  the  El  function  key’s  display  to  "Restart",  and  presents  the  following  message  to  the 
user:  "Processing  Paused".  To  continue  the  processing  and  return  the  screen  to  its  initial 
state,  the  user  need  only  press  the  El  function  key  a  second  time.  The  F5-STOP  function  key 
can  be  pushed  at  any  time  and  presents  the  user  with  the  following  dialog  window:  "Stop 
processing  (Y/N)?".  The  default  response  is  No  and  returns  the  user  to  the  previous  state  of 
the  program,  while  a  Yes  response,  stops  all  processing,  closes  all  files,  and  returns  the  user 
to  the  Wave  Selection  window. 


As  the  METDATA  file’s  records  are  processed,  the  results  are  saved  to  three  ASCII 
files:  RFCOND  or  IRCOND,  TDPxxx,  MP_H_xxx,  and  CN2_xxx.  The  RFCOND  or 
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TABLE  Vn 


General  format  for 

the  RFCOND  or  IRCOND  outuut  files 

Line# 

Variables 

Comment/Example 

1 

5  integers  (RFCOND) 

5  int.  and  1  real  (IRCOND) 

-day,  month,  hour,  minute,  measurement  interval  (min.),  and 
wavelength  (micron) 

-wavelength  is  only  given  for  the  IRCOND  file 

20  2  8  20  5  -  RFCOND  example 

20  2  8  20  5  0.55  -  IRCOND  example 

The  following  two  lines  are  output  for  each  weather  set. 

2 

1 1  real  numbers  for  each 
weather  condition 

-air  temperature  measurement  height  (m),  atmospheric  pressure 
(mbar),  air  tenperature  (“  C),  relative  humidity  (%),  solar 
radiation,  original  wind  speed  (before  any  step  changes),  wind 
direction,  water  temperature  (“  C),  wind  speed  measurement 
height  (m),  wave  height  (m),  and  visibility  (km). 

-  example 

15.0  1010.0  15.0  80.0  500.0  5.0  270  10.0  20.0  1.5  9.0 

2 

9  real  numbers  for  each  weather 
condition 

-stability  error  code,  duct  height  (m),  Monin-Obukhov  length 
(m),  the  three  scaling  parameters  u.  (m/s),  0„.  (“C)  and  q, 

(kg/kg),  and  the  three  roughness  lengths  (m),  Zoe,  (m)  and 
z^(m). _ 

IRCOND  file  is  generated  depending  upon  whether  the  calculation  is  for  the  radar  domain 
or  the  infrared  (visible)  domain,  respectively.  The  format  for  either  of  these  two  files  is 
given  in  Table  Vn.  The  stability/error  code  provides  the  user  with  an  indicator  as  to  the 
reliability  of  the  results.  The  format  of  this  code  is  'IJ.IC,  where  J  can  have  a  value  between 
0  and  8  depending  upon  how  many  of  the  entered  parameters  are  tested  (not  all  parameters 
are  tested)  and  found  to  be  outside  their  respective  minimum  and  maximum  bounds  (see 
Table  VIII).  Thus  J  equal  to  5  would  indicate  that  five  of  the  parameters  were  found  to  be 
outside  their  bounds.  Those  found  to  be  outside  their  bounds  are  set  to  their  respective 
default  values.  I  equal  to  2  indicates  that  a  solution  was  found  by  incrementing  the  wind 
speed  in  steps  of  0.5  m/s,  I  equal  to  1  that  no  solution  could  be  found  up  to  a  maximum  wind 
speed  of  25  m/s,  and  I  equal  to  0  that  a  solution  was  found  without  having  to  change  the  wind 
speed.  If  a  solution  is  found,  K  equal  to  1  indicates  that  the  weather  condition  was  stable, 
and  a  value  of  2  that  the  weather  condition  was  unstable. 
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TABLE  Vm 


Parameter  Bounds  and 

Defaults 

Parameter/Units 

Minimum  Bound 

Maximum  Bound 

Default 

Wind  Speed  (m/s) 

0 

not  tested 

1 

Sea  Temperature  (“C) 

-10 

40 

15 

Air  Temperature  (“C) 

-40 

40 

15 

Wave  Height  (m) 

0 

5 

0.5 

Wind  Speed  Measurement  Height  (m) 

0 

not  tested 

10 

Temperature  Measurement  Height  (m) 

0 

not  tested 

10 

Air  Pressure  (mb) 

0 

not  tested 

1010 

Relative  Humidity  (%) 

0 

100 

80 

The  files  TDPxxx,  MP_H_xxx,  and  CN2_xxx  are  produced  for  each  and  every 
weather  condition  where  xxx  is  an  index  referring  to  the  xxxth  weather  condition  in  the 
input  METDATA  file.  The  TDPxxx  files  contain  vertical  profiles  for  the  air  temperature 
(°C),  relative  humidity  (%),  total  air  pressure  (mbar),  and  wind  speed  (m/s),  the  MP_H_xxx 
files  contain  vertical  profiles  of  the  modified  reffactivity,  and  its  first  and  second  derivatives, 
and  the  CN2_xxx  files  contain  vertical  profiles  for  the  refractive  index,  temperature,  water 
vapour  density,  and  temperature-water  vapour  density  structure  parameters.  All  the  profiles 
go  up  to  max_height  as  given  by  Table  H.  The  modified  refractivity  (M)  is  defined  by 


M(z)  =  Z(z)  +  0.157Z 

(51a) 

MXz)  =  Z'{z)  +0.157 

(51b) 

M"{z)  =  Z"{z) 

(51c) 

where  z  is  the  height  in  meters  and  Z  is  the  refractivity.  Table  DC  gives  the  general  format 
used  by  these  three  files. 

3.2  BATCH  Programs 


If  either  of  the  two  BATCH  programs  are  compiled,  the  user  can  run  the  resultant 
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TABLE  IX 

General  format  for  the  TDPxxx.  MP  H  xxx.  or  CN2  xxx  output  files 


Line# 

Variables 

Comment/Example 

1 

Information  lines  as 
column  titles 

-  MP_H_xxx  files 
#h[m],M(h),M’(h),  M”(h) 

-  TDPxxx  files 

#  HEIGHT  TEMPERATURE  Rel .  Hum.  PRESSURE  WIND  SPEED 

#  [m]  [Deg]  [%]  [mbar]  [m/s] 

#  - -  ■  .  .  . . . . 

-  CN2_xxx  files 

#  HEIGHT  INDEX  TEMP.  VAP.  DEN.  TEMP.-VPDN 

#  [m]  [m-2/3]  [m-2/3]  [m-2/3]  [m-2/3] 

#  . — - -  ===  ■  . — - 

2... 

real  numbers  for  each 
weather  condition 

-  MP_H_xxx  files 

height,  modified  refractivity,  1st  and  2nd  derivatives 

-  TDPxxx  files 

height,  air  temperature,  relative  humidity,  total  air  pressure  and  wind 
speed 

-  CN2_xxx  files 

height,  refractive  index,  temperature,  water  vapour  density  and 
temperature-water  vapour  density  structure  parameters 

programs  as  DOS  command  line  functions.  As  mentioned  earlier,  when  it  is  compiled  with 
the  WKDCN2  definition  parameter  set  to  0,  the  program’s  principal  objective  is  to  produce 
output  files  containing  vertical  profiles  of  the  refractive  index  and  its  gradients.  Likewise, 
when  it  is  compiled  with  the  WKDCN2  definition  parameter  set  to  1 ,  the  program’ s  principal 
objective  is  to  produce  output  files  containing  vertical  profiles  of  the  structure  parameters. 
In  fact,  the  programs  obtained  by  both  of  these  compilation  modes  follow  essentially  the 
same  program  structure  as  the  INTERACTIVE  program  (see  Fig.  3a)  when  it  is  used  in 
BATCH  mode.  Apart  from  the  fact  that  the  input  file  must  be  named  “METDATA”,  and  that 
the  file  CN2_xxx  is  not  produced  when  the  WKDCN2  definition  parameter  is  set  to  0,  and 
that  only  the  CN2_xxx  file  is  produced  when  the  WKDCN2  definition  parameter  is  set  to  1, 
the  differences  are  slight.  The  remaining  differences  have  to  do  with  the  display  of  Error  and 
Status  information,  and  whether  or  not  this  information  is  sent  to  the  screen,  or  to  one  of  the 
two  output  files  “MsgList.txt”  or  “LogFile.txt”.  These  operational  differences  are  controlled 
by  five  programming  switches  (M,  A,  L,  X,  &  C)  that  can  be  added  to  the  command  line. 
The  M  switch  puts  the  program  into  a  non-verbose  mode  (i.e.,  fewer  messages  are  sent  to  the 
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screen  or  to  the  files).  The  A  switch  tells  the  program  to  append  all  Error  &  Status 
information  to  any  previously  existing  output  files  instead  of  creating  new  output  files.  The 
L  switch  directs  the  program  to  redirect  all  screen  output  to  the  file,  “LogFile.txt”,  and  to 
open  the  file  if  necessary  (dependent  upon  the  A  switch).  If  the  program  stops  prematurely, 
the  X  switch  directs  the  program  to  output  any  relevant  messages  to  the  screen,  regardless 
of  the  other  switches.  Finally,  the  C  switch  directs  the  program  to  check  if  a  previous  run 
of  the  program  exited  prematurely.  It  does  this  by  checking  the  previous  “MsgList.txt”  file; 
consequently,  it  requires  the  A  switch  in  order  to  function.  If  it  finds  that  it  did  exit 
prematurely,  the  program  stops  immediately. 
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INTERACTIVE 


BATCH 


FIGURE  3a  -  Schematic  of  the  L(W)WKD  program’s  operating  structure 
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FIGURE  3b  -  Schematic  of  the  L(W)WKD  program’s  operating  structure 
after  selection  of  the  MODIFY  menu 


FIGURE  3c  -  Schematic  diagram  of  the  L(W)WKD  program’s  operating  structure  after 
selection  of  the  Output  Data  menu 
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[f- . =  Weather  Measurements  ====^ 

1  MEASURED  PARAMETERS  -  1 

jf====  Weather  Profiles  =====] 

PROFILE  PARAMETERS 

Horizontal  Range  (km) 

0.000  km 

Wind  Rough  Height 

2.680E-05  m 

Mean  Wave  Height  (Hl/3) 

0.250  m 

Wind  Scale 

9.815E-02  m/s 

Wind  Sensor  Height 

15.000  m 

Wind  Speed 

5.000  m/s 

Mean  Wind  Speed 

5.000  m/s 

Temp.  Rough  Height 

3.837E-05  m 

T/P/RH  Sensor  Height 

15.000  m 

Temperature  Scale 

1.055E-01  C 

Mean  Air  Temperature 

10.000  C 

Air  Temperature 

10.000  C 

Sea  Surface  Temperature 

5.000  C 

Spec.  Hum.  Rough  Ht 

3.066E-04  m 

Mean  Atmospheric  Pressure 

1010.000  mb 

Spec.  Humidity  Scale 

1.756E-05 

Mean  Relative  Humidity 

0.800 

Relative  Humidity 

0.800 

Sea  Surface  Rel.  Humidity 

0.982 

Vapour  Pressure 

9.814  mb 

Vapour  Pressure 

9.814  mb 

Monin-Obukhov  Length 

6.557  m 

Scaling  Pot.  Tang?.  Pressure 

Number  of  Iterations 

8 

KONDO  Tot.  1000  mb 

Fl-  Interactive  P2-  Batch  PS-  ****  F4-  ****  P5-  EXIT 

FIGURE  4  -  L(W)WKD’s  Main  screen  with  theProgram  menu 


PI-  Modify  P2-  Load  Data  P3-  Save  Data  P4-  Output  Data  P5-  RETURN 


FIGURE  5  -  L(W)WKD’s  "Other  Weather  Measurements"  window  with  the  Main  menu 


PI-  Up  P2-  Down  P3-  Edit  P4-  ***  P5-  STOP  P6  -  End 


FIGURE  6  -  The  "Warnings",  "Status"  and  "METDATA  File"  windows  used  in  "Batch" 
mode  with  the  Main  menu 
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4.0  EXAMPLES 

The  objectives  of  this  chapter  are  to  illustrate  some  of  L(W)WKD’s  capabilities  when 
used  as  an  interactive  program.  Examples  of  some  of  the  graphics  that  can  be  created  using 
L(W)WKD  in  interactive  mode  are  shown. 

Figures  7  to  10  show  various  graphs  of  the  type  of  data  that  can  be  produced  by 
version  7.09  of  the  L(W)WKD  program  (non-wavy  model)  for  visible  light  (say  550  nm). 
For  each  graph,  the  water  temperature  is  taken  to  be  12.5  °C,  the  wind  speed  is  taken  to  be 
5  m/s  at  a  height  of  22.5  m  above  the  water  level,  the  relative  humidity  is  taken  to  be  65%, 
and  the  atmospheric  pressure  to  be  1026  mbar  at  a  height  of  12  m  above  the  water  level.  The 
Hi;3  wave  height  is  taken  to  be  0. 18  m,  and  the  air  temperature  measured  at  12  m  above  the 
water  level  is  taken  to  be  8.4  °C  when  the  ASTD  is  -4. 1  (unstable  case),  and  16.6  when 
the  ASTD  is  4.1  °C  (stable  case). 

Figure  7  shows  a  plot  of  the  vertical  air  temperature  (symbols  only)  and  refractivity 
(solid  line  with  symbols)  profiles  produced  under  both  of  these  conditions.  Comparing  the 
vertical  air  temperature  profiles  for  the  cases,  one  notices  that  the  gradient  near  the  water 
surface  (below  1  m)  is  much  stronger  in  the  unstable  case  than  for  the  stable  case  and  of 
opposite  sign.  But  at  higher  elevations,  the  gradient  is  stronger  in  the  stable  case  than  in  the 
unstable  case  and  of  the  same  sign.  Comparing  the  vertical  refractivity  profiles,  one  notices 
the  same  type  of  behaviour,  except  that  there  is  actually  an  inversion  in  the  refractivity  for 
the  unstable  case.  This  implies  that,  for  the  stable  case,  light  rays  are  continuously  bent 
towards  the  surface  of  the  earth  with  the  strongest  bending  occurring  closest  to  the  surface; 
while  in  the  strongly  unstable  regime,  light  rays  are  bent  slightly  towards  the  surface  of  the 
earth  above  ~6  m  and  away  from  the  surface  of  the  earth  below  ~6  m.  These  are  the  kinds 
of  conditions  that  can  cause  ducting  of  the  light,  and  secondary  image  formation, 
respectively. 
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Air  Temperature  &  Refractivity  Vertical  Profiles 

Water  Temp,  m  12.5  °C;  RH  >  65%  «  12  m;  Wind  Speed  >  5m/e  O  22.6  m 


I  A  ASTD» -4.10:  Temp,  o  ASTD  -  4.1C;  Temp. -*-ASTD  . -4.10;  Ret. -e-ASTD  .  4.1C;  Ref.  | 
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FIGURE  7  -  Plots  of  the  vertical  air  temperature  and  visible  refractivity  profiles  for 
ASTDs  of  -4. 1  °C  and  4. 1  °C  using  the  L(W)WKD  model  in  non-wavy  mode 


Refractivity  Structure  Parameter  Vertical  Profile 

Water  Temp.  >  12.5  °C;  RH  «  65  %  e  12  m;  Wind  Speed  ^  5m/s  O  22.6  m 


Height  (m) 

I-  ■  —  •  ASTD  -  .4.10;  NW  — ASTD  -  4.1C:  NW  | 

FIGURE  8  -  Plot  of  the  refractivity  structure  parameter  for  ASTDs  of  -4. 1  °C  and  4. 1  °C 
using  the  L(W)WKD  model  in  non-wavy  mode 
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structure  Parameter  vs  ASTD  @  2  m  above  MWL 

Water  Temp.  =  12.5  °C;  RH  =  65  %  6  12  m;  Wind  Speed  =  5m/s  S  22.6  m 


Refractivitv . Temperature  — I — Vapour  Density  ~  .X- -Temp.-Vap.  Density  | 


FIGURE  9  -  Plot  of  the  refractivity  structure  parameter  and  its  three  components  as  a 
function  of  changing  air-sea  temperature  difference  (ASTD)  using  the 
L(W)WKD  model  in  non-wavy  mode 


Virtual  Potential  Temperature  Profile  Parameters  vs  ASTD 

Water  Temp.  =  12.5  ®C:  RH  =  65  %  ©  12  m;  Wind  Speed  =  5m/s  ©  22.6  m 


Scale  Param.  -  H — Roughness  Length  \ 


FIGURE  10  -  Plot  of  the  scale  parameter  and  roughness  height  for  the  virtual  potential 
temperature  as  a  function  of  air-sea  temperature  difference  (ASTD)  using  the 
L(W)WKD  model 
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Figure  8  shows  a  plot  of  the  calculated  vertical  refractivity  structure  parameter 
profiles  for  these  same  unstable  and  stable  cases.  As  can  be  seen,  the  structure  parameter  for 
the  unstable  case  is  predicted  to  be  as  much  as  5  times  greater  than  that  for  the  stable  case 
at  heights  below  ~4.5  m,  and  to  become  increasingly  weaker  than  the  stable  case  as  the 
height  becomes  increasingly  greater  than  ~4.5  m.  In  fact,  at  a  height  of  ~38  m  it  is  10  times 
less.  The  reason  for  this  is  due  to  the  different  forms  of  the  structure  parameter  functions  that 
are  used  in  the  unstable  and  stable  regimes  (see  Eq.  46). 

Figure  9  shows  a  plot  of  the  calculated  refractivity  structure  parameter,  and  its  three 
components,  as  a  function  of  the  air-sea  temperature  difference  at  a  height  of  2  meters  above 
the  mean  water  level  (MWL)  for  the  same  conditions.  As  can  be  seen,  the  temperature 
structure  parameter  has  the  form  of  a  parabola  with  a  minimum  near  an  ASTD  of  zero 
degrees,  and  the  vapour  density  stmcture  parameter  decreases  in  a  linear  manner  from  about 
0.2  for  an  ASTD  of -10  °C  to  a  minimum  close  to  an  ASTD  of  6  °C  before  increasing  slightly 
to  0.03  at  an  ASTD  of  10  °C.  The  temperature-vapour  density  structure  parameter  decreases 
almost  linearly  from  a  value  of  0. 18  at  an  ASTD  of  -10  °C  to  a  value  of  zero  near  an  ASTD 
of  0.5  °C,  increases  slightly  to  about  0.02  at  an  ASTD  of  4  °C,  before  decreasing  to  -0.08  for 
an  ASTD  of  10  °C.  The  appropriate  summation  of  these  three  structure  parameters  results 
in  the  refractivity  structure  parameter  (solid  line).  Looking  closely  at  the  curve,  one  notices 
that  it  is  always  greater  than  the  air  temperature  structure  parameter  for  negative  ASTD,  and 
less  than  the  air  temperature  strucmre  parameter  for  ASTDs  greater  then  6  °C. 

Figure  10  shows  a  plot  of  the  calculated  virtual  potential  temperature  scaling 
parameter  (solid  line)  and  roughness  height  (dashed  line)  as  a  function  of  the  ASTD  for  these 
same  conditions.  The  scaling  parameter  is  seen  to  increase  linearly  from  about  -0.4  for  an 
ASTD  of -10  °C  to  zero  near  an  ASTD  of  0  ®C,  and  then  to  increase  almost  linearly  to  about 
0. 15  for  an  ASTD  of  10  ”C.  On  the  other  hand,  the  roughness  height  has  a  shape  reminiscent 
of  the  arctangent  or  Heaviside  step  function,  with  values  from  23  to  26  pm  for  ASTDs  below 
-1.5  °C,  and  from  35  to  38  pm  for  ASTDs  above  1.5  °C. 
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5.0  MAPTIP  Comparisons 

This  section  shows  new  calculations  of  the  maximum  inter- vision  range  (MIVR)  and 
minimum  mirage  range  (MMR)  using  L(W)WKD  as  a  module  within  IRBLEM  compared 
with  earlier  results  and  the  measurements  made  during  the  1993  MAPTIP  trial  (Ref.  2).  The 
Marine  Aerosol  Properties  and  Thermal  Imager  Performance  (MAPTIP)  trial  was  held  in  The 
Netherlands  from  Oct.  18  to  Nov.  3, 1993,  and  was  sponsored  by  RSG  8  of  NATO’s  panel  4. 
The  principle  goals  of  the  trial  were  to  study  the  concentrations  and  properties  of  marine 
aerosols,  and  to  study  the  performance  of  imaging  systems  within  the  marine  boundary  layer 
(MBL)  (Ref.  1 5).  One  of  the  studies  (see  Ref.  2),  involved  measuring  the  minimum  mirage 
ranges  (MMRs)  and  maximum  inter-vision  ranges  (MIVRs)  for  a  set  of  halogen  lamps 
installed  at  various  heights  aboard  the  Hr.  Ms.  Tydeman  as  it  sailed  away  from  a  set  of 
visible  and  infrared  cameras  placed  at  the  beach  station  in  the  town  of  Katwijk.  Under 
unstable  conditions,  the  MIVR,  for  an  observer  (camera)  and  a  source  (lamp)  at  given 
elevations,  is  the  range  at  which  the  source  goes  beyond  the  observer’s  horizon,  while  the 
MMR  is  the  range  at  which  a  secondary  image  (or  mirage)  of  the  source  first  appears. 

Tables  XI  to  XIV  show  all  the  measured  MMR  and  MTVR  MAPTIP  data  for  cases 
A  thru  M  obtained  using  our  Visible  1  and  Visible  2  camera  systems.  They  also  contain 
prediction  results  obtained  from  calculations  using  Version  7.09  of  the  L(W)WKD  model, 
and  the  earlier  version  discussed  in  Ref.  1  (Ver.  6.0).  The  colunms  in  each  of  these  tables, 
starting  with  the  first  column,  show  the  measurement  case  (A  to  M),  the  name  of  the  light 
source,  the  height  of  the  source  above  the  water  level  and  its  uncertainty  (discussed  later), 
the  experimentally  measured  MMRs  or  MIVRs  and  a  measure  of  their  uncertainties,  the 
earlier  LWKD  (non-wavy  mode)  predicted  MMRs  or  MIVRs  and  the  difference  between 
them  and  the  measured  values,  the  earlier  LWWKD  (with  the  wavy  mode)  predicted  MMRs 
or  MTVRs  and  the  difference  between  them  and  the  measured  values,  and  the  LWKD  Ver. 
7.09  (non-wavy  mode)  predicted  MMRs  or  MIVRs  and  the  difference  between  them  and  the 
measured  values.  The  average  and  standard  deviation  for  each  of  the  difference  colunms 
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(model  prediction  -  measurement)  is  given  at  the  bottom  of  each  table  along  with  the  number 
of  measurements  used  in  the  calculations.  Looking  at  these  values,  it  appears  that,  on 
average,  the  new  MMR  predictions  are  at  least  as  good  as  those  obtained  using  the  earlier 
LWWKD  predictions,  and  that  the  new  MTVR  predictions  are  about  0.3  km  less  than  both 
the  earlier  LWKD  and  LWWKD  predictions. 

At  this  point,  I  would  like  to  come  back  to  the  estimation  of  the  uncertainty  in  the 
height  of  each  of  the  light  sources.  Looking  at  Tables  XI  to  XTV,  one  notices  that  these 
values  vary  with  respect  to  both  the  particular  measurement  Case,  and  the  particular  Object 
( light  source).  This  is  due  to  the  pitching  motion  of  the  Hr.  Ms.  Tydeman  (like  a  teeter- 
totter)  that  is  generated  by  the  waves.  Assuming  that  the  amplitude  of  the  pitching  motion 
for  each  of  the  light  sources  is  only  dependent  upon  the  wave  height,  I  have  defined  the 


Normalized  Elevations  (Source  Height/Hi/3)  for  Each  Light  Source 


PTGURE 1 1  -  Plot  of  the  normalized  elevations  for  each  light  source  and  for  MAPTIP  cases 
C,  F,  H,  K,  M,  and  their  average 
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following  parameter,  which  I  call  the  normalized  elevation  For  each  source  n,  it  is 

defined  as: 


H 


f/"  - 

•^max  ^min 


H 


(52) 


1/3 


where  H1/3  is  the  H1/3  wave  height,  and  and  are  the  minimum  and  maximum 
measured  heights  of  each  source  above  the  ships  water  line.  From  our  video  recordings,  for 
several  of  the  MAPTIP  cases,  it  was  possible  to  make  a  number  of  measurements  of 
and  Figure  1 1  shows  the  results  for  each  of  the  visible  lights  obtained  from  an  analysis 
of  MAPTIP  cases  C,  F,  H,  K,  and  M.  The  average  over  all  the  cases  is  given  by  the  solid  line 
and  the  value  of  each  average  is  located  next  to  each  point.  Despite  the  large  scatter  that  can 
exist  between  the  data  for  different  cases,  it  is  clear  that  the  normalized  elevations  are  larger 
for  the  lower  elevation  lights  and  much  smaller  for  the  higher  lights.  The  reason  for  this  is 
due  to  the  positioning  of  the  lights  on  the  ship.  The  lowest  lights  are  positioned  at  the  ship’s 
stem,  while  the  highest  light  is  on  its  mast  near  the  axis  about  which  the  ship  pitches,  and  the 
other  lights  are  somewhere  in  between.  Thus,  using  the  average  values  for  the  normalized 
elevations,  the  uncertainties  in  the  source  heights  given  in  Tables  XI  to  XTV  for  each  source, 
6„,  are  obtained  using  the  following  relationship; 

^  H,„V1  ,  (53) 


where  <Hn„onn>  is  the  mean  value  for  each  source  given  in  Table  X.  Thus  we  can  define  a 
lower  and  upper  limit  for  each  source  height  given  by; 


where  Hn„,^ 
(no  waves). 


r/"  _  JJ’'  _  S 

^lower  ^mean 

H"  =  ht’'  +8 

^  upper  ^mean 


(54) 


is  the  measured  height  of  each  source  above  the  ships  water  line  for  a  flat  sea 
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TABLE  X 


Average  Normalized  Elevations  for  each  Light  Source 


Light 

IL/IR 

2C/2Ca 

3I73R 

Chimney 

4C 

1.43 

1.16 

0.66 

0.24 

0.24 

Tables  XV  to  XVin  show  the  results  obtained  by  running  L(W)WKD  Version  7.09 
in  the  non-wavy  mode  for  three  different  source  heights  (mean,  lower  limit,  and  upper  limit), 
and  two  different  values  for  von  Karman’s  constant  for  each  of  the  MAPTff  cases.  The 
mean  source  height  for  each  source  is  that  given  in  Tables  XI  to  XFV,  and  the  lower  and 
upper  limits  are  calculated  using  Eq.  54.  For  each  table,  column  one  gives  the  MAPTIP  case, 
column  two  the  light  source,  and  column  three  the  measured  MMR  or  MIVR.  The  following 
sets  of  columns  show  the  results  of  the  model  calculations  for  von  Karman  constants  of  0.40 
and  0.41  for  the  mean,  lower  limit  and  upper  limit  source  heights.  In  addition.  Tables  XVI 
and  XVin  (contain  MIVRs),  also  give  the  results  for  a  non-refractive  atmosphere  (i.e.  for 
free-space  propagation).  The  calculation  using  a  value  of  0.41  for  von  Karman’s  constant 
(default  value  is  0.40)  were  made  due  to  the  results  recently  reported  by  Dion  and  Gardenal 
(Ref.  26).  They  found  that  they  obtained  a  better  fit  for  radio  frequency  (RF)  data  to  model 
calculations  using  refractivity  profiles  generated  by  L(W)WKD  in  non-wavy  mode  when  a 
value  of  0.41  was  used  instead  of  0.40. 

Figures  12  to  15  are  graphical  representation  of  the  data  in  Tables  XVI  to  XVin. 
Figure  12  shows  the  MMR  results  obtained  using  Ver.  7.09  of  the  model  in  non-wavy  mode 
(x),  and  Ver.  6.0  in  wavy  mode  (o)  for  source  heights  at  their  mean  heights  above  the  water 
level  versus  the  experimental  data.  The  dashed  line  and  the  solid  line  give  their  respective 
least  squares  fits  and  their  respective  equations  are  given  at  the  top  of  the  graph.  While  there 
is  a  scatter  of  about  1  km  about  the  straight  lines,  and  the  fits  are  quite  good  for  both  models, 
our  newer  model  does  perform  slightly  better.  The  MMRs  predicted  by  the  new  model  are 
about  1.5%  less  than  the  measured  values  while  those  predicted  by  the  earlier  model  are 
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Model  Estimates  vs  Measurements  for  Minimum  Mirage  Ranges 
MAPUP  Visible  1  &  2  Cameras 


X  LWKDV7.09-Mean  o  LWWKD  V6.0  -  Mean 

. Unear  (LWKDV7.09- Mean)  - Unear  (LWWKD  V6.0  •  Mean) 

FIGURE  12  -  Plot  of  L(W)WKD  Ver.  7.09  results  using  the  non-wavy  mode  (x),  and  Ver. 

6.0  using  the  wavy  mode  (o)  versus  experimental  data  obtained  from  our 
Visible  1  and  Visible  2  cameras  during  the  MAPTIP  trial.  The  solid  and 
dashed  lines  give  the  respective  least  squares  fits.  The  mean  height  for  the 
light  sources  are  used  in  each  case 

about  2.8%  greater. 

Figure  1 3  shows  the  MMR  results  obtained  using  Ver.  7.09  of  the  model  in  non-wavy 
mode  with  a  von  Karman  constant  of  0.40  (x),  and  with  a  von  Karman  constant  of  0.41  (o) 
for  source  heights  at  their  lower  limits  above  the  water  level  versus  the  experimental  data. 
The  dashed  line  and  the  solid  line  give  their  respective  least  squares  fits  and  their  respective 
equations  are  given  at  the  top  of  the  graph.  Again  there  is  a  scatter  of  about  1  km  about  the 
straight  lines,  and  both  fits  are  very  good  and  definitely  better  than  the  model  calculations 
given  in  Fig.  12.  Looking  at  the  fits,  one  notices  that  the  MMRs  predicted  using  both  values 
of  von  Karman ’s  constant  would  precisely  predict  the  measured  MMRs  if  it  were  not  for 
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Model  Estimates  vs  Measurements  for  Minimum  Mirage  Ranges 
MAPTIP  Visibie  1  &  2  Cameras 


X  LWKDV7.09- Lower  {K=0.40)  o  LWKDV7.09- Lower  (K=0.41) 

- Unear  (LWKD  V7.09  •  Lower  (K=0.40)) . Unear  (LWKD  V7.09  -  Lower  (K=0.41)) 

FIGURE  13  -  Plot  of  L(W)WKD  Ver.  7.09  results  using  the  non-wavy  mode  for  von 
Karman  constants  equal  to  k  =  0.40  (x)  and  k  =  0.41  (o)  versus  experimental 
data  obtained  from  our  Visible  1  and  Visible  2  cameras  during  the  MAPTIP 
trial.  The  solid  and  dashed  lines  give  the  respective  least  squares  fits.  The 
lower  limits  for  the  heights  of  each  light  source  is  used  in  each  case 

constant  offsets  of  0.44  km  and  0.26  km,  respectively.  In  other  words,  for  a  von  Karman 

constant  of  0.40,  the  model  consistently  predicts  an  MMR  that  is  0.44  km  shorter  than  the 

measured  value,  while  for  a  constant  of  0.41  it  consistently  predicts  a  value  that  is  0.26  km 

shorter  than  that  measured.  Consequently,  a  value  of  0.41  for  von  Karman’ s  constant  along 

with  the  use  of  the  lower  limits  for  the  source  heights  gives  the  best  fit  to  the  measured  MMR 

data. 


Figure  14  shows  the  MIVR  results  obtained  using  Ver.  7.09  of  the  model  in  non-wavy 
mode  (x),  Ver.  6.0  in  wavy  mode  (o),  and  for  a  non-refractive  atmosphere  (□)  for  source 
heights  at  their  mean  heights  above  the  water  level  versus  the  experimental  data.  The  dashed 
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Model  Estimates  vs  Measurements  for  Maximum  Intervision  Ranges 
MAPT1P  Visibie  1  &  2  Cameras 


X  LWKDV7.09-Mean  o  LWWKD  V6.0  -  Mean  ■  Non-Ref. 

. Linear  (LWKD  V7.09  -  Mean)  - Linear  (LWWKD  V6.0  -  Mean) - Linear  (Non-Ref.) _ 

FIGURE  14  -  Plot  of  L(W)WKD  Ver.  7.09  results  using  the  non-wavy  mode(x),  Ver.  6.0 
using  the  wavy  mode  (o),  and  free  space  calculations  (□)  versus  experimental 
data  obtained  from  our  Visible  1  and  Visible  2  cameras  during  the  MAPTIP 
trial.  The  solid  and  dashed  lines  give  the  respective  least  squares  fits.  The 
mean  height  for  the  light  sources  are  used  in  each  case 

line  and  the  solid  line  give  the  respective  least  squares  fits  for  the  two  L(W)WKD  model 

calculations  and  their  respective  equations  are  given  at  the  top  of  the  graph.  Again  there  is 

a  scatter  of  about  1  km  about  the  straight  lines,  and  while  the  fits  are  quite  good  for  both 

L(W)WKD  models,  our  newer  model  does  perform  slightly  better,  and  both  perform 

substantially  better  than  a  prediction  using  a  non-refractive  atmosphere.  The  MTVRs 

predicted  by  the  new  model  both  have  offsets  of  approximately  1.5  km.  Thus  both  models 

tend  to  underestimate  the  measured  MTVR  data  by  up  to  1.5  km,  with  the  difference 

diminishing  for  larger  MIVRs. 


Figure  15  shows  the  MTVR  results  obtained  using  Ver.  7.09  of  the  model  in  non-wavy 
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mode  with  a  von  Karman  constant  of  0.40  (x),  and  with  a  von  Karman  constant  of  0.41  (o) 
for  source  heights  at  their  lower  limits  above  the  water  level  versus  the  experimental  data. 
The  dashed  line  and  the  solid  line  give  their  respective  least  squares  fits  and  their  respective 
equations  are  given  at  the  top  of  the  graph.  Again  there  is  a  scatter  of  about  1  km  about  the 
straight  lines,  and  both  fits  are  very  good  and  definitely  better  than  the  model  calculations 
given  in  Fig.  14.  Looking  at  the  fits,  one  notices  that  the  MTVRs  predicted  using  both  values 
of  von  Karman’s  constant  would  precisely  predict  the  measured  MIVRs  if  it  were  not  for 
constant  offsets  of  0.63  km  and  0.50  km,  respectively.  In  other  words,  for  a  von  Karman 
constant  of  0.40,  the  model  consistently  predicts  an  MIVR  that  is  0.63  km  shorter  than  the 


Model  Estimates  vs  Measurements  for  Maximum  Intervision  Ranges 
MAPTIP  Visibie  1  &  2  Cameras 


X  LWKDV7.09- Upper  (K=0.40)  o  LWKD  V7.09  -  Upper  (K=0.41) 

. Unear  (LWKD  V7.09  -  Upper  (K=0.40)) - Unear  (LWKD  V7.09  -  L^jper  (K=0.41)) 

FIGURE  15  -  Plot  of  L(W)WKD  Ver.  7.09  results  using  non-wavy  mode  for  von  Karman 
constants  equal  to  k  =  0.40  (x)  and  k  =  0.41  (o)  versus  experimental  data 
obtained  from  our  Visible  1  and  Visible  2  cameras  during  the  MAPTIP  trial. 
The  solid  and  dashed  lines  give  the  respective  least  squares  fits.  The  upper 
limits  for  the  heights  of  each  light  source  is  used  in  each  case 
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measured  value,  while  for  a  constant  of  0.41  it  consistently  predicts  a  value  that  is  0.50  km 
shorter  than  that  measured.  Consequently,  a  value  of  0.41  for  von  Karman’s  constant  along 
with  the  use  of  the  lower  limits  for  the  source  heights  gives  the  best  fit  to  the  measured  MMR 
data. 


Thus,  for  the  models  studied  in  this  report.  Figs.  12  to  15  show  that  the  best  fit  to  the 
measured  MMR  and  MIVR  MAPTIP  data  is  obtained  for  Ver.  7.09  of  the  L(W)WKD  model 
in  non-wavy  mode,  when  the  source  heights  are  corrected  for  the  pitching  of  the  ship,  and 
von  Karmen’s  constant  is  taken  to  be  0.41.  Breaking  this  down,  we  see  that  correcting  for 
the  ships  pitching  motion  resulted  in  fits  of  the  MMR  and  MIVR  experimental  data  to  the 
model  results  that  are  almost  perfectly  linear  such  that  only  slight  y-offsets  of  approximately 
0.5  km  remained.  This  implies  that  the  model,  on  average,  only  underestimates  the  measured 
values  by  approximately  0.5  km.  Consequent  modification  of  von  Karman’ s  constant  to  0.4 1 
from  0.40  produced  a  decrease  in  the  y-offset  by  about  0.2  km,  implying  that  increasing  the 
values  to  0.42  or  0.43  could  lead  to  an  even  better  fit. 

6.0  CONCLUSIONS 

This  report  has  described,  in  detail,  both  the  physics,  and  mathematical  structure 
contained  within  the  most  recent  version  (Version  7.09)  of  the  L(W)WKD  model,  and  how 
it  compares  with  the  previously  published  version  (Ref.  1),  and  the  PIRAM  model  (Ref.  3). 
It  has  also  given  a  complete  description  of  the  operational  procedures  required  to  operate  the 
associated  L(W)WKD  computer  program.  It  also  presented  and  discussed  a  number  of 
results  which  can  be  obtained  from  the  model. 

A  comparison  between  earlier  measurements  taken  during  the  MAPTIP  trial  and  the 
earlier  L(W)WKD  model  calculations  shown  in  Ref.  2  with  more  recent  calculations  using 
Version  7.09  of  the  L(W)WKD  has  also  been  performed.  This  new  analysis  also  involved 
a  modification  to  the  source  heights  to  account  for  the  pitching  of  the  ship,  and  a  modifica- 
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tion  to  von  Karman’s  constant  from  0.40  to  0.41.  With  these  two  modifications,  the  least 
square  fits  of  the  MMR  and  MIVR  experimental  data  to  the  model  results  that,  on  average, 
only  underestimates  the  measured  values  by  less  than  0.5  km. 

Future  improvements  to  the  model,  and  to  the  program  include  the  option  of  using 
other  functional  forms  for  the  gradient  functions,  and  structure  parameter  functions.  This 
could  also  involve  the  possibility  of  using  different  values  for  the  various  gradient 
coefficients,  and  other  parameters  that  are  currently  fixed  by  extending  the  number  of 
variables  defined  by  the  initialization  file  “LWWKD.IM”.  Additionally,  it  would  be  useful 
to  integrate  the  calculation  of  aerosol  vertical  profiles,  and  to  allow  files  containing  this 
profile  to  be  created  by  the  program  when  it  is  being  operated  in  BATCH  mode.  The 
production  of  more  meteorological  profiles  and  structure  profiles  may  also  be  useful. 
Finally,  further  study  and  perhaps  modification  of  the  refractivity  profile  (perhaps  by 
modifying  von  Karman’s  constant)  may  be  required  in  order  to  better  fit  the  MMR  and 
MIVR  predictions  given  using  Ver.  7.09  of  L(W)WKD  model  with  measurements  obtained 
during  the  MAPTBP  trial. 
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TABLE  XI 

Calculated  and  Measured  Minimum  Mirage  Ranges  tMAPTIP  -  Visible  1  Camera 
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Calculated 


TABLE  XI  (continued) 

and  Measured  Minimum  Mirage  Ranees  CMAPTIP  -  Visible  1  Cameral 


CaseG  CaseF  CaseE  CaseD  CaaeC  i  CaseB 
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TABLE  Xn 

Calculated  and  Measured  Maximum  Intervision  Ranges  ('MAPTIP  -  Visible  1  Camera 
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TABLE  Xn  (continued') 

Calculated  and  Measured  Maximum  Intervision  Ranees  rMAPTIP  -  Visible  1  Camera) 


Case  F  Case  E  Case  D  Case  C  Case  B  Case  A 
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TABLE  Xm 

Calculated  and  Measured  Minimum  Mirase  Ranees  fMAPTIP  -  Visible  2  Camera) 
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TABLE  XDI  (continued) 

Calculated  and  Measured  Minimum  Mirage  Ranees  CMAPTIP  -  Visible  2  Camera) 
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TABLE  XIV 

Calculated  and  Measured  Maximum  Intervision  Ranges  fMAPTIP  -  Visible  2  Camera 
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TABLE  XIV  (continued) 

Calculated  and  Measured  Maximum  Intervision  Ranees  fMAPTP  -  Visible  2  Camera') 
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TABLE  XVn 

Calculated  and  Measured  MMRs  for  Different  Source  Heights  IMAPTIP  -  Visible  2 
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TABLE  XVm 

Calculated  and  Measured  MTVRs  for  Different  Source  Heieht  fMAFTP  -  Visible  21 


UNCLASSIFIED 

84 


7.0  ACKNOWLEDGEMENTS 

I  would  like  to  thank  the  efforts  of  Jacques  Claverie  of  France  and  Yvonnick  Hurtaud 
of  CELAR  in  France  for  their  many  useful  discussions  while  we  worked  towards  a  common 
vision,  and  even  more  concretely,  towards  common  physical  equations  for  use  in  current  and 
future  maritime  boundary  layer  models.  I  would  also  like  to  thank  Denis  Dion  of  DREV  and 
Lionel  Gardenal  of  Informission  for  their  help  and  advice  during  the  various  stages  of  this 
project. 


8.0  REFERENCES 


1.  Forand,  J.L.,  “The  L(W)WKD  Marine  Boundary  Layer  Model”,  DREV-R-9618, 
March  1997,  UNCLASSIFIED 

2.  Forand,  J.L.,  Dion,  D.,  Hurtaud,  Y.,  and  Stein,  K.,  “MAPTIP  Work  Group  Report: 
Refractive  Effects  in  the  Visible  and  the  Infrared”,  DREV  R-9621,  June  1997, 
UNCLASSIFIED 

3.  Claverie,  J.,  Hurtaud,  Y.,  de  Fromont,  Y.  and  Junchat,  A.,  "Modelisation  des  Profils 
Verticaux  d'Indice  de  Refraction  et  de  C\  en  Atmosphere  Marine",  Propagation 
Assessment  in  Coastal  Environments,  AGARD  Conference  Proceedings  567,  p.  29- 
1,  1994. 

4.  Monin,  A.S.  and  Obukhov,  A.M.,  "Basic  Regularity  in  Turbulent  Mixing  in  the 
Surface  Layer  of  the  Atmosphere",  Trad.  Geophys.  Inst.  ANSSSR,  No.  24,  p.  163, 
1954. 

5.  Lumley,  J.L.  and  Panofsky,  H.A.,  "The  Structure  of  Atmospheric  Turbulence", 
Interscience  Monographs  and  Texts  in  Physics  and  Astronomy,  Vol.  Xn,  p.  99, 
Edited  by  R.E.  Marshak,  University  of  Rochester,  Rochester,  New  York,  1964. 

6.  Hill,  R.J.,  Clifford,  S.F.  and  Lawrence,  R.S.,  "Refractive-Index  and  Absorption 
Fluctuations  in  the  Infrared  caused  by  Temperature,  Humidity  and  Pressure 
Fluctuations",  J.  Opt.  Soc.  Am.,  Vol.  70,  No.  10,  pp.  1192-1205, 1980. 

7.  Hill,  R.J.  and  Lawrence,  R.S.,  "Refractive  Index  of  Water  Vapor  in  Infrared 
Windows",  Infrared  Phys.,  Vol.  26,  No.  6,  pp.  371-376, 1986. 


UNCLASSMED 

85 


8.  Dion,  D.  and  Schwering,  P.,  “On  the  Analysis  of  Atmospheric  Effects  on  Electro- 
Optical  Sensors  in  the  Marine  Surface  Layer”,  IRIS  Conference,  UK,  June  1996 

9.  Low,  T.B.  and  Hudak,  D.R.,  “Final  Report  on  the  Development  and  Testing  of  a 
Marine  Boundary  Layer  Model”,  KelResearch  Corp.  Report  under  DSS  contract 
#W7701-8-2419/01-XSK,  September  1990 

10.  Beaulieu,  A.J.,  "Atmospheric  Refraction  Model  and  the  Effects  of  Surface  Waves", 
DREV  R-4661/92,  May  1992,  UNCLASSIFIED 

1 1 .  Walmsley,  J.L.,  "On  Theoretical  Wind  Speed  and  Temperature  Profiles  over  the  Sea 
with  Applications  to  Data  from  Sable  Island,  Nova  Scotia",  Atmosphere  Ocean,  Vol. 
26,  pp.  203-233,  1988. 

12.  Kondo,  J.,  "Air-Sea  Bulk  Transfer  Coefficients  in  Diabetic  Conditions",  Bound. 
Layer  Meteo.,  9, 1975,  pp  91-1 12. 

13.  Claverie,  J.,  "Private  communication",  July  1998. 

14.  Businger,  J.A.,  Wyngaard,  J.C.,  Isumi,  Y.  and  Bradley,  E.F.,  "Flux-Profile 
Relationships  in  the  Atmospheric  Surface  Layer",  J.  Atmos.  Sci.,  28,  pp.  181-189, 
March  1971. 

15.  Smith,  S.D.,  "Water  Vapour  Flux  at  the  Sea  Surface",  Boundary-Layer  Meteorology, 
47,  pp.  277-293, 1989. 

16.  Garratt,  J.R.,  "The  Atmospheric  Boundary  Layer",  Cambridge  University  Press, 
Cambridge,  p.  286, 1992. 

17.  Bom,  M.,  and  Wolf,  E.,  "Principles  of  Optics",  6th  (Corrected)  Ed.,  Pergamon  Press, 
Toronto,  pp.  87-98, 1980. 

18.  Edlen,  B.,  "The  Dispersion  of  Standard  Air",  J.  Opt.  Soc.  Am.,  43,  pp.  339-344, 
1953. 

19.  Erickson,  K.E.,  "Investigation  of  the  Invariance  of  Atmospheric  Dispersion  with  a 
Long-Path  Refractometer",  J.  Opt.  Soc.  Am.,  52,  pp.  777-780, 1962. 

20.  Bean,  B.R.  and  Dutton,  E.J.,  "Radio  Meteorology",  National  Bureau  of  Standards 
Monograph  92,  p.  1-22, 1966. 

21.  Edson,  J.B.,  Fairall,  C.W.,  Mestayer,  P.G.  and  Larsen,  S.E.,  "A  Study  of  the  Inertial- 
Dissipation  Method  for  Computing  Air-Sea  Fluxes",  J.  of  Geoph.  Res.,  96,  pp. 


UNCLASSMED 

86 


10689-10711, 1991. 

22.  Wyngaard,  J.C.,  Izumi,  Y.  and  Stuart  Jr.,  A.S.,  "Behavior  of  the 
Refractive-Index-Stmcture  Parameter  Near  the  Ground",  J.  Opt.  Soc.  Am.,  61,  pp. 
1646-1650, 1971. 

23.  Bataille,  P.,  "Analyse  du  comportement  d'un  systeme  de  telecommunication  optique 
fonctionnant  a  0.83  pm  dans  la  basse  troposphere",  these  de  doctorat  de  TUniversite 
de  Rennes,  Nov.  1, 1992. 

24.  Davidson,  K.L.,  Schacher,  G.E.,  Fairall,  C.W.,  and  Goroch,  A.K.,  "Verification  of 
the  Bulk  Method  of  Calculating  Overwater  Optical  Turbulence",  Applied  Optics,  20, 
pp.  2919-2924, 1981. 

25.  Forand,  J.L.,  "Horizontal  Non-homogeneous  Effects  in  the  Marine  Boundary  Layer", 
Propagation  and  Imaging  through  the  Atmosphere,  SPIE  Meeting,  Denver  CO,  pp. 
180-191, 1997. 

26.  Dion,  D.  and  Gardenal,  L.,  "On  the  accuracy  of  bulk  methods  for  predicting  path 
losses  in  the  marine  surface  layer"  ,1999  IEEE  AP-S/URSI  International  Symposium, 
Orlando,  USA,  July  11-16, 1999. 

27.  Hurtaud,  Y.,  Dion,  D.,  Claverie,  J.  and  Forand,  L.,  "Cooperations  Franco-canadienne: 
Arrangement  Specifiques  n°7  sur  la  propagation  eiectromagnetique  au-dessus  de  la 
mer.  Rapport  Final",  DGA  GEOS/99.49/SOP/80832,  le  13  avril  1999. 


UNCLASSIFffiD 

87 


APPENDIX  A 

PHYSICAL  CONSTANTS  AND  VARIABLES 


Ouantitv 

Symbol  Value 

Units 

Height  above  mean  water  level  (m) 

z 

m 

Wind  speed  measurement  height 

Zo 

m 

Temp.,  pressure,  rel.  humidity  measurement  height 

Zl 

m 

Wind  speed 

u 

m/s 

Wind  speed  constant  of  integration 

Uo 

m/s 

Wind  speed  roughness  length 

Zou 

m 

Wind  speed  scaling  constant 

u. 

m/s 

Wind  speed  gradient  function 

9u 

Integrated  wind  speed  gradient  function 

'T 

^  U 

Sea  temperature 

T 

K 

Virtual  sea  temperature 

Tvs 

K 

Virtual  potential  sea  temperature 

0VS 

K 

Air  temperature 

T 

K 

Mean  virtual  temperature  in  MBL 

T 

vm 

K 

Virtual  temperature 

Tv 

K 

Virtual  potential  temperature 

ev 

K 

Mean  virtual  potential  temperature  in  MBL 

Qvn, 

K 

Virt.  pot.  air-sea  temperature  difference  (VPASTD) 

A0v(z,) 

K 

Virt.  pot.  temperature  constant  of  integration 

0VO 

K 

Virt.  potential  temperature  roughness  length 

m 

Virtual  potential  temperature  scaling  constant 

Bv. 

K 

Virtual  potential  temperature  gradient  function 

V 

Integrated  virt.  pot.  temperature  gradient  function 

¥ 

Specific  humidity 

q 

g/g 

Air-sea  specific  humidity  diference  (ASSH) 

Aq(Zi) 

g/g 
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Specific  humidity  constant  of  integration 
Specific  humidity  roughness  length 
Specific  humidity  scaling  constant 
Specific  humidity  gradient  function 
Integrated  specific  humidity  gradient  function 
Acceleration  of  gravity 
Von  Karman’s  constant 
Monin-Obukhov  length 
Normalized  height  (=  z/L) 

WKD  stable  gradient  function  coefficient  ^  0) 
WKD  stable  wind  speed  gradient  function  coefficient 
WKD  stable  virt.  pot.  temp,  gradient  function  coef. 
WKD  stable  specific  humidity  gradient  function  coef. 
Unstable  gradient  function  coefficient  <  0) 

Unstable  wind  speed  gradient  function  coefficient 
Unstable  virt.  pot.  temp,  gradient  function  coef. 
Unstable  specific  humidity  gradient  function  coef. 
Kondo  gradient  function  coefficient  ^  0) 

Kondo  stable  wind  speed  gradient  function  coefficient 
Kondo  stable  virt.  pot.  temp,  gradient  function  coef. 
Kondo  stable  specific  humidity  gradient  function  coef. 
Molecular  mass  of  dry  air 
Molecular  mass  of  water 

Ratio  of  mol.  masses  of  water  and  dry  air  (=  M^/MJ 

Boltzmann's  constant 

Universal  gas  constant 

Dry  air  gas  constant  (=  R/M^) 

Water  vapour  gas  constant  (=R,yM,^,) 


fio 

^Oq 

q* 

q 

g  9.81 

K  0.4 

L 


C 

a 

a, 

02 

P 

P. 

P2 

P3 

y 

Yi 

72 

73 

M, 

M 

k 

R 

K 

Rw 


5 

5 

5 

5 
16 
16 
16 
16 

6 
6 
6 
6 

28.9 

18.0 

0.622 

1.381x10-2^ 

8.31 

0.288 

0.462 


g/g 

m 

g/g 

m/s2 

m 


g 

g 

J/K 

J/mol/K 

J/g/K 

J/g/K 
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Moist  air  gas  constant 

Rn. 

J/g/K 

Total  atmospheric  pressure 

P 

mbar 

Total  atmospheric  reference  pressure 

Pref 

1000 

mbar 

Dry  air  atmospheric  pressure 

Pa 

mbar 

Water  vapour  pressure 

Pw 

mbar 

Total  air  density 

P 

g/m^ 

Dry  air  density 

Pa 

g/m' 

Water  vapour  density 

Pw 

g/m' 

Mixing  ratio 

W 

g/g 

Spec,  heat  capacity  of  dry  air  at  constant  pressure 

s 

1.004 

J/g/K 

Spec,  heat  capacity  of  dry  air  at  constant  volume 

Cv 

0.717 

J/g/K 

Ratio  of  heat  capacities  for  dry  air  (=cjc^) 

Y 

1.40 

Saturated  water  vapour  pressure  over  water 

P  ' 

mbar 

Saturated  water  vapour  pressure  over  salt  water 

Psw' 

mbar 

Relative  humidity 

H. 

Dew-point  temperature 

Ta 

K 

Salinity 

S 

g/kg 

Salinity  coefficient 

8 

5.37x10'^ 

10-m  neutral  heat  flux  coefficient 

CtN 

1.0x10'^ 

10-m  neutral  moisture  flux  coefficient 

^EN 

1.2x10'^ 

Chamock’s  constant 

«c 

0.011 

Dynamic  viscosity 

mVs 

‘a^’  coefficient  in  wind  dependent  relation 

9.267x10-8 

m^/s/^C 

‘bj’  coefficient  in  wind  dependent  relation 

K 

1.346x10-5 

mVs 

‘C’  coefficient  in  wind  dependent  relation 

c 

-5.5 

Average  of  the  1/3  highest  wave  heights 

H,;3 

m 

Water  wave  height 

Hw 

m 

Water  wavelength 

Xw 

m 
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Molecular  refractivity  for  dry  air 

Da 

Molecular  refractivity  for  water  vapour 

Dw 

Refractivity  coefficient  for  dry-air 

A 

K/mbar 

Refractivity  coefficient  for  water  vapour 

B 

K/mbar 

i*  Line  strength 

Si 

Wavelength 

X 

pm 

i*  wavelength  resonance 

pm 

Refractive  index 

n 

Refractivity  [=  (n  - 1)  x  10"^] 

Z 

Refractivity  for  dry  air 

Za 

Refractivity  for  water  vapour 

Zw 

Refractive  index  structure  parameter 

m-2^ 

Refractivity  index  structure  parameter 

m'^ 

Temperature  structure  parameter 

Cji 

m-2^' 

Temp.-specific  humidity  structure  parameter 

^Tq 

m-^3 

Specific  humidity  structure  parameter 

m-2/3 

Virtual  temperature  structure  parameter 

Cj  2 

m-2/3 

Virtual  temp.-spec.  humidity  structure  parameter 

Virt.  potential  temperature  structure  parameter 

Virt.  pot.  temp.-specific  humidity  structure  parameter 

-2/3 

m 

Virt.  pot.  temp.-water  vap.  density  structure  parameter 

-2/3 

Water  vapour  density  structure  parameter 

-2/3 

m 

i*  structure  parameter  function 

fi 

Structure  parameter  function  ratio  (=f2/fi) 

A21 

1 

Virt.  pot.  temp.-spec.  humidity  correlation  coefficient 

±0.8 
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APPENDIX  B 

LISTING  OF  MODIFICATIONS 


//  Mar.  1995 

//  Ver.  5.0  Oldest  version. 
//  Sept.  1995 


//  Ver.  6.0 
// 

// 

//  June  1996 
//  Ver.  6.01 
// 

//  Nov.  1996 
//  Ver.  6.11 
// 

// 

// 

// 

// 

// 

// 

// 

// 

//  Nov.  1996 
//  Ver.  6.2 
// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


1)  Modified  to  include  batch  file  option  (Sept.  95)  by 
Lionel  Gardenal  -  Informission  Inc. 

2)  Added  copyright  screen. 

1)  Modified  to  use  the  hydrostatic  equation  to  calculate  the  total 
pressure  instead  of  the  dry-air  pressure. 

1)  use  dry-air  or  total  air  pressure  for  potential  temp,  calculation 

2)  enter  desired  adiabatic  pressure  for  potential  temp,  calculation 

3)  can  now  choose  between  LWWKD  or  Kondo  profile  for  stable  case 
The  defaults  are  LWWKD,  Dry-Air  and  1000  mb. 

Note:  while  these  options  can  not  be  directly  chosen  in  the  “batch" 
mode,  if  the  user  goes  into  the  "interactive"  mode,  modifies  them  and 
then  returns  and  runs  the  “batch"  mode,  the  batch  mode  will  use  these 
modified  options  until  they  are  again  changed  in  "interactive"  mode. 

4)  the  saved  binary  file  is  slightly  different  as  it  saves  the  options; 
thus,  the  file  ID  is  now  1. 

Modified  to  include  multiple  weather  parameter  sets  where  each  one  is 
distinguished  by  a  horizontal  range. 

1)  the  horizontal  range  is  a  new  parameter  which  has  been  added  to  the 
PARAM  Structure. 

2)  there  are  two  new  options  (Add  &  Delete)  in  the  modify  weather  menu. 

The  'Add'  )cey  adds  a  new  weather  set  which  can  subsequently  be  edited, 
and  the  'Delete'  key  removes  the  currently  displayed  weather  set. 

3)  the  user  can  add  up  to  10  weather  sets  and  can  cycle  between  them  using 
the  LEFT  and  RIGHT  arrows.  These  keys  can  also  be  used  in  other  windows. 
For  example,  if  you  are  in  the  graphics  menu,  change  the  weather  set  which 
is  displayed  to  get  graphs  pertaining  to  it. 

4)  the  weather  sets  are  always  arranged  in  order  of  increasing  range. 

5)  the  saved  binary  file  is  again  different  (file  ID  =  2)  as  it  saves 
information  for  all  the  weather  sets.  The  file  saves  all  the  weather 
set  data  as  well  as  the  refraction  data  for  a  chosen  wavelength  and 
model  as  before,  except  that  now  it  has  it  for  all  the  weather  sets. 

This  means  that  the  file  can  now  have  range-dependent  refractivity  data. 

6)  Note:  the  options  are  the  same  for  all  weather  sets  and  the  wave  height 
for  all  the  sets  are  always  the  same.  This  makes  the  file  easier  to  use 
in  ray  tracing. 

7)  Added  extra  plot  so  that  different  profile  parameters  can  be  plotted 
against  the  adiabatic  pressure  (900  -  1200  mb) . 


//  April  17,  1997 

//  Ver  6.21  1)  Corrected  problem  with  structure  parameter  calculations.  The  calculation 

//  of  pow(x,2/3),  etc.  were  giving  1.  Now  changed  to  pow(x, 2 . 0/3 . 0) . 


//  May 
//  Ver 
// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


20,1997 

6.22  1)  Corrected  problem  in  the  subroutines  wvap_struct  &  vptemp_wvap_struct 

of  the  module  profile.c.  The  problem  was  that  the  temperature  was  being 
converted  to  Kelvin  when  it  shouldn't  have  been.  This  caused  there  to  be 
two  conversions  from  Celsius  to  Kelvin;  one  in  the  above  functions  and 
a  second  in  the  subroutines  vapour_density  and  air_density  of  the  module 
weather. c.  These  corrections  only  affect  the  structure  parameter  calcu¬ 
lations. 

2)  A  second  problem  was  noticed  with  respect  to  the  calculation  of  the 
virtual  potential  temperature  structure  parameter,  specific  humidity 
structure  parameter,  and  the  virtual  potential  temperature-specific  humidity 
structure  parameter.  A  factor  of  K^2  (von  Karmen's  constant)  must  be  added 
to  each  of  the  subroutines  virt_pot_temp_struct,  spec_humid_struct,  and 
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// 

// 

// 

// 

// 

// 

// 

/ /  Dec . 
//  Ver. 
// 

// 

//  Apr. 
//  Ver. 
// 

// 

// 


vpteinp_shuin_struct  of  the  module  profile.c.  This  is  due  to  an  error  in  my 
LWWKD  report. 

3)  Noticed  a  problem  with  the  calculate_weather (mode)  subroutine  in  the 
graphics  modules.  The  problem  was  that  they  were  not  using  the  proper 
value  for  mode.  The  calls  to  calculate_weather  were  changed  from 
calculate_weather (mode)  to  calculate_weather (pot_mode) .  "pot_mode"  is  a 
global  constant. 

19,  1997 

6.3  1)  Allows  the  possibility  that  the  surface  relative  humidity  can  be 

something  other  than  100% 

2)  Corrected  error  in  the  Kondo  formula  (5.0  to  6.0) 

02,  1998 

6.31  1)  Correct  calculation  of  duct  height  in  batch  section  of  the 

program . 

2)  Correct  problem  where  an  IRCOND  file  was  created  for  the  radar 
case  instead  of  the  desired  file  RFCOND 


//  July 
//  Ver. 
// 

// 

// 

/ /  Sept . 
//  Ver. 
// 

// 

// 

// 

n 

II 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


16, 

6.40 


28, 

7.00 


1998 

1)  Modify  certain  modules  so  that  most  input/output  functions 
are  separate  from  any  calculation  functions.  This  will  ma)ce  it 
easier  to  incorporate  the  core  LWWKD  routines  into  other 
operating  systems . 

1998 

1)  Modifications  are  made  to  a  number  of  routines  in  order  to 
comply  with  common  options  (options  commune)  agreed  upon  t>y 
Canada  &  France. 

a)  The  vertical  pressure  profile  is  calculated  using  the  same 
earlier  equation  except  that  the  virtual  temperature  is  used. 

b)  In  the  equation  that  relates  the  wind  speed  roughness  length 
to  the  wind  speed  scaling  constant,  the  v  term  (dynamic  viscosity) 
is  now  a  function  of  the  virtual  potential  temperature. 

c)  The  relative  humidity  at  the  surface  has  a  default  value  of 
98.2%.  This  is  good  for  salt  water  with  a  salinity  of  34g/)cg 
(RHsea  =  1  -  5.37e-3  x  S) . 

d)  The  Kondo  function  is  now  the  default  mode.  This  means  that 
the  Kondo  function  is  used  for  the  stable  case  (VPASTD  >  0) . 

e)  The  Monin-Obu)diov  length  is  now  calculated  using  the  average 
of  the  virtual  temperature  instead  of  the  average  temperature. 

f)  The  scaling  parameters  have  been  modified  so  that  they  agree 
with  the  definitions  used  by  the  French.  In  other  words,  our 
scaling  equations  have  been  multiplied  by  von  Karman's  to  give 
those  used  by  the  French. 

2)  In  order  to  ma)ce  these  changes,  a  number  of  changes  were  also 
made  to  many  of  the  subroutines,  including  changing  their  names, 
moving  them  to  a  different  or  newly  created  module,  and  rewriting 
their  codes . 

3)  Made  slight  changes  to  ASCII  &  Postscript  files  containing  the 
Met  &  Profile  data. 

4)  The  program  can  be  compiled  in  either  batch  or  interactive 

by  setting  the  INTERACT  definition  to  (0)  or  (1)  respectively.  In 
batch  mode  the  program  can  easily  be  interfaced  with  any  other 
set  of  programs  that  are  able  to  provide  it  with  a  'metdata'  file. 

5)  The  batch  part  of  the  program,  the  LWWKD. C  &  LWWKDBAT.C  modules, 
owes  many  of  its  modifications  to  the  wor]c  of  Lionel  Gardenal  of 
Le  Groupe  Inf omission  Inc.  The  two  modules  have  been  changed  so 
that  they  produce  a  file  containing  the  vertical  refractive  index 
profiles  (MP_H_xxx)  for  wavelengths  in  visible,  infrared  or  radar 
bands  under  given  weather  meteorological  conditions  and  sea  state. 
It  also  produces  a  file  containing  vertical  profiles  of  main 
meteorological  parameters  (TDPxxx) .  Finally,  it  creates  a  file 
(RFCOND  or  IRCOND)  containing  values  of  the  internal  parameters  of 
the  used  model.  This  latter  is  the  Wavy  WKD  boundary  layer  model 
developed  at  DREV.  Evaluation  of  atmospheric  stability  is  based 
the  value  of  the  Monin-Obu)chov  length. 

a)  The  following  run-time  arguments  have  been  added  to  LWWKD. C 
-A  :  Append;  The  log  file  infomation  will  also  be  appended 
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// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


if  in  Log  mode  (see  -L) .  If  A  is  not  specified, 

1. the  previous  output  messages  file  will  be  erased. 

2.  Log  file  will  be  first  reseted. 

-C  :  Conditional;  If  in  Append  mode  and  if  a  previous  program 
has  already  been  exited  this  one  will  also  be  stopped. 

No  effect  if  not  in  Append  mode  since  the  message  file 
will  first  be  erased. 

-L  :  Log;  Information  normally  sent  to  the  output  screen  will 
be  sent  to  a  log  file  instead.  No  interaction  with  user. 
Abnormal  execution  messages  will  be  sent  to  the  output 
message  file  and  log  file,  (see  also  -A  arg.) 

b)  Modifications  of  1997-04-03  (Lionel  Gardenal,  Informission) 
Profile  of  water  vapour  density  is  replaced  by  the  profile  of 

relative  humidity  in  TDPxxx  output  file. 

c)  Modifications  du  1997-07-11  (Lionel  Gardenal,  Informission) 

When  the  input  values  are  outside  valid  range,  LWWKD  displays 
a  Monin-Obukhov  length  equal  to  9999.99;  furthermore,  the 
TDPxxx  file  contains  constant  values  for  the  meteorological 
profiles.  These  values  are  equals  to  measured  data  as  read 
by  LWWKD  from  the  METDATA  file. 

d)  Modifications  du  1997-07-22  :  (Lionel  Gardenal,  Informission) 
TDPxxx  files  contains  scalar  values  of  measured  wind  speed 
(instantaneous  and  mean  value  on  24h) ,  visibility  and  air 
mass  parameter.  These  parameters  are  used  by  MERGPROF. 

e)  Modifications  du  1997-07-30  :  (Jean  Gill>ert,  Informission) 
Adding  of  running  options  to  the  program  permetting  to  redirect 
output  to  a  log  file  and  to  send  copy  abnormal  messages  to  a 
message  file. 

f)  Modifications  du  1998-08-19  :  (Lionel  Gardenal,  Informission) 

A  new  process  has  taeen  integrated.  It  is  called  when  the  WKD 
model  doesn't  converge.  It  consists  of  increasing  the  wind 
speed  until  a  solution  is  obtained  (convergence) . 


//  Sept.  29,  1998 

//  Ver.  7.01  1)  Added  a  new  graph  type  so  that  one  can  now  plot  the  dependence 

//  of  the  fit  parameters  against  changes  in  von  Karmen's  constant. 

//  2)  Added  a  new  parameter  to  the  PARAM  structure.  The  new  parameter 

//  is  the  wind  speed  at  the  surface.  For  the  moment  it  is  set  to  zero 

//  and  can  not  be  changed  from  within  the  program.  The  only  way  to 

//  change  it  is  by  recompiling  the  program. 

//  Oct.  9,  1998 

//  Ver.  7.02  1)  Changed  B  coefficient  for  radar  to  373256/T. 

//  2)  Corrected  A  coefficients  used  with  the  refractivity  structure  parameter. 

//  Previous  values  were  77.7,  77.5,  77.5  &  78.5  instead  of  78.5,  77.7,  77.5  & 

//  77.5  for  the  VIS,  NIR,  MIR  &  FIR  respectively. 

//  3)  Modified  calculation  of  duct  height  used  in  Ijatch  mode.  Added 

//  a  new  routine  "duct_height"  and  modified  ’'save_mp_h" . 

//  Oct.  16,  1998 

//  Ver.  7.03  1)  Corrected  stability/error  output  in  batch  mode  when  a  fit  is  found 

//  by  modifying  the  wind  speed. 

//  Oct.  27,  1998 

//  Ver.  7.04  1)  Modified  the  RF  refractivity' s  B  coefficient  to  agree  with  Bean  &  Dutton 

//  Thus  the  B  coefficient  is  now  '72  +  3.75E5/T'. 

//  Nov.  13,  1998 

//  Ver.  7.05  1)  Corrected  error  in  loading  a  previous  created  "Iwwjtd"  saved  file. 

//  Mar.  3,  1999 

//  Ver.  7.06  1)  Corrected  error  noticed  by  Espen  Lippert,  of  Norway,  pertaining  to  our 

//  B  coefficients.  As  a  result  we  have  changed  the  values  in  this  version. 

//  The  change  is  less  than  1%  in  the  visible,  near  and  mid  IR  and  about  20% 

//  in  the  far  IR. 

//  May  4,  1999 

//  Ver.  7.07  1)  Add  output  of  a  file  containing  structure  parameter  profiles  (C2N,  etc.) 

//  2)  Change  structure  parameter  graphs  so  that  they  allow  the  user  to  choose 

//  the  x-axis  range  and  the  number  of  data  points.  */ 

//  3)  Corrected  values  of  coindex  output  by  "save_mp_h"  when  their  is  no  fit. 

//  4)  A  second  compiling  option  has  been  added  to  the  program  for  operation  in 
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// 

// 

// 

// 

// 

// 

// 

//  June 
//  Ver. 
// 

// 

//  June 
//  Ver. 
// 

// 

// 

// 

// 


batch  mode.  By  setting  the  WKDCN2  definition  to  (1),  the  program  is 
compiled  so  that  it  functions  like  the  WKDCN2  routine  in  IRBLEM.  If  set  to 
(0)  it  is  compiled  so  that  it  functions  like  the  LWWKD  routine  in  IRBLEM. 
Note:  this  compiling  option  is  only  valid  when  it's  compiled  in  BATCH  mode. 
Note:  can  not  use  the  wavy  model. 

5)  For  INTERACTIVE  compilation  the  BATCH  mode  also  produces  profiles  of  the 
the  structure  parameters. 

10,  1999 

7.08  1)  Corrected  error  in  Cn2  calculation.  The  program  was  subtracting  the 

the  cross-term  Ctq  instead  of  adding  it.  I.E.  the  -2  was  changed  to  +2. 

2)  Also  changed  all  C++  style  comments  to  C  style  comments. 

15,  1999 

7.09  1)  Mike  Duffy  added  a  user-configurable  parameter  file  (LWWKD.INI).  It 

allows  the  modification  of  Kel_cel,  Mair,  Mwater,  R,  g,  gamm,  cp,  Kconst, 
p_adiab,  max_height,  pot_mode,  and  scale_mode.  More  will  be  added  at  a 
later  date. 

2)  The  routine  Runmode  was  moved  from  LWWKDBAT.C  to  LIBERROR.C.  Other  minor 
changes  were  also  made  to  Mike's  original  code. 
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